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ABSTRACT 


Main purpose of this study is to determine the accuracy of the 
line spring model, to investigate the effect of interaction between 
two and three cracks, and to provide extensive numerical results 
which may be useful in applications. Line spring model with 
Reissner’3 plate theory is formulated to be used for any number and 
configurations of cracks provided that there is symmetry. This model 
is used to find stress intensity factors for elliptic internal 
cracks, elliptic edge cracks and two opposite elliptic edge cracks. 
Unfortunately, because of the unavailability of previous work done 
on the cases considered, only stress intensity factors for central 
elliptic crack could be compared with other methods. Despite the 
simplicity of the line spring model, the results are found to be 
surprisingly close. 





1. INTRODUCTION 


From the viewpoint of practical applications, the analysis of a 
part through crack in a structural canponent which may locally be 
represented by a "plate" or a "shell" is certainly one of the most 
important problems in fracture mechanics. In its general form, the 
problem is a three-dimensional crack problem in a bounded ge one try 
where the stress fields perturbed by the crack interacts very 
strongly with the surfaces of the solid. At present, even for the 
linearly elastic solids, a neat analytical treatment of the problem 
very heavily rely on some kind of numerical technique, such as 
alternating method, [5], [6], or boundary integral equation method, 
[7]; but most notably on the finite element method, [8]-[10]. The 
renewed interest in reoent years in the so called "line-spring 
model" first described in [3] has been due partly to the desire of 
providing simpler and less expensive solutions to the part-through 
crack problem and partly to the fact that for certain important 
crack geometries, the model seems to give results that have an 
acceptable degree of accuracy. 

In a plate or shell containing a part- thro ugh crack and 
subjected to membrane and bending loads, the net ligament(s) around 
the crack would generally have a constraining effect on the crack 
surface displacements and rotations. The basic idea underlying the 
"line-spring model" consists of approximating the three-dimensional 
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crack problem by a two-dimensional coupled bending-membrane problem 
through the reduction of the net ligament stresses to the neutral 
surface of the plate or shell as a membrane load N and a bending 
moment M. In the resulting two-dimensional problem, the crack 
surface displacements are represented by a crack opening 
displacement 6 and a crack surface rotation ®, referred to, again, 
the neutral surface. The quantities N, M, $ and 8 are assumed to be 
functions of a single variable, namely the coordinate x-j , along the 
crack in the neutral surface (Fig. 1). The pair of functions (&,&) 
or (N,M) are determined from the corresponding mixed boundary value 

problem for the plate or the shell having a through crack in which N 

\ 

and M are treated as unknown crack surface loads. Once N and M are 
determined, the stress intensity factors are evaluated from the two- 
dimensional elasticity solution of a strip under the membrane force 
N and the bending moment M. 

The model introduced in [3] is based on the classical theory. 
However, the asymptotic stress field around the crack tip given by 
classical plate bending theory is not oonsistant with the elasticity 
solutions, whereas a transverse shear theory (such as that of 
Reissner's) which can accomodate all stress and moment resultants on 
the crack surface seperately ( i. e. , three boundary conditions in 
plates, five in shells) give results which are identical to the 
asymptotic solutions obtained from the plane strain and anti-plane 
shear crack problems [2], [11]. The line spring model was later used 
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in [12], [13] to treat the longitudinal part- through crack problem 
in a cylinder by using the classical shell theory and in [1], [14] 
by using transverse shear theory. Rather extensive results for 
corner cracks and for collinear surface cracks in a plate having a 
finite width are given in [15]. The similar problem of surface 
crack-boundary interaction in a cylindrical shell with a free or a 
fixed boundary is considered in [16] and [17]. 

In this study, stress intensity factors for elliptic internal 
cracks under pure bending and tension, and for collinear elliptic 
edge cracks under pure tension will be studied by using the line- 
spring model with Reissner's plate bending theory. Interaction of 
two and three identical elliptic edge and internal creaks will also 
be .studied by using the same method. 
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2. FORMULATION OF THE PROBLEM 


2. 1 Governing Integral Equations 

The problem under consideration is a surface or an internal 
crack problem for a relatively thin-walled structural component 
which is solved basically as a plate or shell problem. In the usual 
notation u.j f u 2 , Ug are the components of the displacement vector, 
$1 and are angles of rotation of the normal to the neutral 
surface in and planes, respectively, and Njj, M^ and 

V^(i, j=1 ,2) are respectively the membrane, moment and transverse 
shear resultants (Fig. 1). 

Related mixed boundary value problem for any nmber of cracks 
has been formulated in [15] by using Reissner's plate bending 
theory. The derivatives of the crack surface displacement and the 
■crack surface rotation on the neutral surface are chosen as the 
unknown functions. Fredholm kernels k^(x,t) and k 2 (x,t) in [15] 
represent the effect of finite length of the plate in x^ direction. 
They will vanish in this case since the plate has an infinite length 
in x 1 direction. Integral equations are written for only half of the 
infinite plate (x^ > 0) and are only valid when symmetry with 

respect to the XgXg plane exists. D represents the normalized 
cracked region on the x 1 axis, which means that if there is a crack 
in the region (b,e), D represents the region (b/a , c/a) where "a" 
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ia the half crack length and equals to (c-b)/2. x is also normalized 
parameter on x^ axis which is xsx^/a^ where a^. the half crack 
length for the corresponding crack(i subscript is needed because 
general formulation is done for more than one crack). 

From [15], integral equations are as follows: 


*(x) 1 f 1 1 

+ | Gl (t) 

L t-x t+x J 


2*w 


dt = 


<r 0 (x) 


, x £ D , (1) 


m(x) 


- Ri 


6E 


1 1 




h [ — * — 1 - % [ 

L t-x t+x J L 


(t^x) 3 (t+x) 3J 


r 1 1 

+ RJ K 2 (p |t-x| ) + Ko(p I t+x| ) 

^ a. 


't-x 


t+x 


G 2 (t) dt 


m Q (x) 

6E 


x £ D , (2) 


where K 2 is the modified Bessel function of the second kind and the 
constants , R 2 , R^, Rj, and p are defined in Appendix I. h is the 
thickness of the plate, 2a is the length and L(x^) is the depth of 
the crack (Fig. 1). E is the modulus of elasticity and p i 3 the 

Poisson's ratio of the material. The unknown functions are defined 
by 
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The external loads 




( 4 ) 


represent uniform membrane and bending resultants applied to the 
plate away from the crack region and <T and m which are defined by 


N(x,0) 

ff(x) 

h 


6M(x,0) 

m(x) = , x£D , 

h 2 


(5) 


are the membrane and bending loads applied to the crack surfaces 
(Fig. 1). The integral equations are obtained from the following 
mixed boundary conditions in X 2=0 plane 



N 22 (x,0) 

= - N 22 * 

N(x) , 

x£D , 

(6a) 

u 2 (x,0) 

= 0 

9 


(6b) 

M2 2 (x,0) 

+ 

CM 

i 

ii 

M(x) , 

, 

(7a) 

* 2 (x,0) 

= 0 

t 

x£D , 

(7b) 


where the general principle of superposition is used to account for 
the loading n| 2 , and m£ 2 applied to the structure away from the 
crack region. 


2.2 Line- Spring Model for an Edge Crack 

Edge crack is chosen for introductory purposes because of the 
simplicity of its formulation. Modifications will be made as more 
complicated crack geometries are examined. 

The configuration studied here is an "infinite" elastic plate 
of thickness h, which contains a surface crack of length 2a and 
depth L(x.j) penetrating part through the thickness (Fig. 1). For 
general purposes L(x 1 )can be any function which enables us to treat 
any crack geometry. At remote distances from the crack site, the 
plate is subjected to loads equipollent to a uniform simple tension 
in the x 2 direction and to pure bending about the x-j axis. The first 
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assumption in the line-spring analysis is that the stress intensity 
factor at a point along the crack front with coordinate x.j is 
identical to the stress intensity factor for an edge cracked strip 
in plane strain (Fig. 2) , subjected to an axial force and moment 
equal, respectively, to N 22 (x 1 ,0} and M 22 (xi,0), and having a crack 
depth equal to L(x^). In [1], [3], mode I stress intensity factor 
for the plane strain problem is given as 


K(s) = )/h [ffg t (s) + mg b (s)] , (8) 

s(x.,) = L(x.j )/h , (9) 


where functions g b (s) and g b (s) are called shape functions for 
tension and bending respectively, and given as 


S t (s) 

Sb< s ) 


= y%S ^b i S 2 ^ i " 1 ^ , 

ini 

r\ 

(10a) 

= /ws Y_ c i s^ i “ 1 ^ 

-L-l 

(10b) 


Ihe value of n is chosen according to the desired accuracy. The 
coefficients b 1 and can be found by a suitable curve fitting 
(Appendix II). Coefficients of the shape functions for a number of 
crack geometries are given in Tables 2, 4, 6, 8, 10, 12, 14, 16, 17. 
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In order to obtain K(x^) and M(x^) in terms of and G 2 , the 
energy available for fracture along the crack front is expressed in 
two different ways, namely as the crack closure energy and as the 
product of load-load point displacement. In a plate with an edge 
crack subjected to uniform tension N and uniform bending moment M 
(Fig. 2), if K is the stress intensity factor given by the plane 
strain solution, from the crack closure energy, the energy (per unit 
width) available for fracture may be obtained as 


73 

G = (U-V) 

c)L 



(ID 


where 0 is the work done by the external loads and V is the strain 
energy. 

_ Let & and © be the load line displacements corresponding to N 
and M shown in Fig. 2 and d& and d© be the changes in & and © as 
the crack length goes from L to L+dL under fixed load conditions. 
Then referring to Fig. 2, the changes in U and V may be expressed as 

dU = NdS + Md© , (12) 

dV = (1/2) [N($+d&) + M(©+d©) - (1/2) (N& + M©) = 

= (1/2) (Nd& + Md© ) , (13) 

Equations (12) and (13) give the energy available for a crack growth 
dL as follows: 
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<1(0- V) = (1/2) (NdS + Md8) . 


(14) 


On the other hand for constant N and M, and for a change of dL in 
the crack length, we have 


a* 

d& = dL , 

3l 


d8 dL . 

3l 


Thus, from (14) and (15) it follows that 


a 

— (0-V) 

Sl 


1 

r 38 Den 

G = — 

N — + M , 

2 

. C)L DL - 


(15) 


(16) 


and, by using (11) we find 


1 c)8 c)e 1-r 2 

— [n — + M — 1 = K 2 

2 L Sl Dl j E 


(17) 


From (5) and (8), we may rewrite (17) in terms of 9 and m as 
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follows: 



as 

C h — + 

a L 


h 2 a s 

6 3l 


1 -9 C 

h (gf«r 

E 


2 s 2 


2g t g b ffm + g b m 2 ) . 


(18) 


In order to solve 9 and m in terms of 8 and • or vice versa we 
introduce the so called compliance coefficients. In a cracked strip, 
displacement and rotation are functions of both the applied stress 
and bending manent. This relation may be expressed as 

hS = + A fcb m , (19a) 

, 1^/6 = A bfc c + A fab m , (19b) 

where (with proper normalizations) A tb =A bt by elastic reciprocity. 
Hie compliance coefficients A^ depend only on L and vanish when 
L=0. If we substitute (19) into (18) and equate common coefficients 
in the quadratic forms of both sides, we get 


dA tt 

dL 


2(1-P 2 )h 

E 




(20a) 
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I 


(20b) 


dA tb ^bt 2(1-l» 2 )b 


dL dL 


•StSb 


dA bb 2(1-l^)h 


dL 


Sb 


By defining 


«ij = h ' 


L 

A J g igj dL , (i, j = b,t) , 


and knowing that A^sO at L=0 , we may write 


2h 2 (1-P 2 ) 

A tt = a tt » 

E 


2h 2 (1-r 2 ) 


A tb = A bt 


a tb » 


2h 2 ( I-? 2 ) 


l bb = 


bb 


(20c) 


( 21 ) 


(22a) 


(22b) 


(22c) 
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Substitution of (22) into (19) gives 


2h(1-P 2 ) 

8 ( a tt ff + c tb“^ » (23) 

E 

I2(i-r 2 ) 

e (a bt <r + cc bb m) , (24) 

E 


where ®tb =a bt’ ^ rom (3)» 8 and 8 may be expressed in terms of the 
unknown functions and G 2 as 


X 

r 


8 = 2 ^ 2 ( x ,+ 0 ) = 2 

G 1 (t) dt , 

(25) 



8 = 2au 2 (x,+0) = 2a 

G 2 (t) dt . 

( 26 ) 


- ! 


If we solve (23) and (24) for <r(x) and m(x) and substitute (25) and 
(26) for 8 and 8, we can determine «r(x) and m(x) in terms of the 
unknown functions G 1 and G 2 as follows: 


/ 


<T(x) = E CY tt (x) G 2 (t) dt + Y fcb (x) 


-I 


/ 


x 

G-j(t) dt] 


(27) 
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X 


m(x) = 6E CY bt (x) 

G 2 (t) 


J 

1 

where 




a 

g bb 

T tt “ 

h(1-F 2 ) 

A 

v - 

1 

g tb 

• tb " 

6(1-? 2 ) 

A 

Y, 

1 

g bt 

T bt 

6h(1-P 2 

r ) A 

ii 

•Q 

XI 

1 

g tt 

36(1-r 2 ) 

A 


G 1 (t) dt] , (28) 


(29a) 


(29b) 


(29c) 


(29d) 


4 ■ - *lb • 


(29e) 


We should keep in mind that o^j is a function of L(x)/h, which means 
is also a function of x. 
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For a single crack shown in Fig. 1 , the normalized cracked 
region on axis (called D in (1) and (2)) is defined as (0,1). (1) 
and (2) then become 


I 


(x) 

(t) dt = — 1. , 0 < x < 1 , ( 30 ) 

E 2*8 ■> L t-x t+x J ' E 




m(x) 


- Ri 


6E j 

o 


1 1 


1 1 


+ - R 3 [ + 1 

t-x t+x *■ (t-x)^ (t+x)3 J 


+ Rj,[ Kg ( p | t-x | ) + Kg(p |t+x| ) j 


t-x 


t+x 


Gg(t ) dt 


m. 


6E 


0 < x < 1 . (31) 


For convenience in the numerical solution of integral equations we 
need to convert the integrals calculated over (0,1) to (-1,1). 
Because of that the geanetry and loading conditions are symmetric 
with respect to the XgX^ plane, Ug(x,+0) and j* 2 (x,+0) should also 
be symmetric with respect to the same plane. In mathemathical form 


u 2 (-x,+0) = u 2 (x,+0) 


(32a) 
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£ 2 (- x * + °) = ^2^ x » +0 ^ 


(32b) 


Observing that derivative of a symmetric function is antisymmetric, 
we may write 


G^(-t) = -G^t) 

t 

(33a) 

4* 

o" 

1 

II 

-U 

1 

• 

(33b) 


By applying the above properties, integrand of (30) can be rewritten 
as 


I 


r 



O 


+ 



dt 



1 

G 1 ( t ) dt 

t-x 



o 


1 

— G.j (t^ ) dt^ 
t^ +x 


(3*0 


By a change of variable t^ = -t and from (33a), we can write 
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t.j +x 


G 1 (t^ ) dt, = 


-t+x 


[-G, (t)] (-dt) 


j 

-I 


t-x 


G, (t ) dt 


( 35 ) 


Substituting (35) in (3*0 


o 


r 1 1 i 

+ 1 G- ( t ) 

t-x t+x J 


dt = 


t-x 


G,(t) dt . (36) 


After applying the same procedure to (30 and substituting (27) and 
(28) for flf(x) and m(x), integral equations take the following final 
form 


Y bb (x) G-j(t) dt - R, 


1 

R 2 


t-x 3 (t-x) 3 


A 

R 4 Kg (p | t-x | ) ] G,(t) dt + *f bt (x) Gg(t ) dt 

t-x J 


m. 


6E 


-I 


, -1 < x < 1 , 


(37a) 
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*tb (x > 


r 

I G-j ( t ) dt + V tt (x) J G 2 (t) 
I -I 

GgCt) 


dt 


1 

2* j 


■ dt = , -1 < x < 1 . 

t-x E 


(37b) 


From (3) it follows that the unknown functions and G 2 must 
satisfy the single valuedness conditions given by 

I I 

J G^t) dt = 0 , J G 2 (t) dt = 0 . (38) 

-I -I 

After solving (37) and (38) by using Gauss-Chebyshev closed type 
integration formula which is defined in Appendix III, a backward 
procedure is applied to find the stress intensity factors along the 
crack front. Unknown functions are integrated (Appendix IV) from -1 
to the value of x at which the stress intensity factor is desired. 
Then those integrals are substituted to (27) and (2 8) with the 
functions Yjj(x), (i, j=t, b) evaluated at that particular point to 
get corresponding <T and m. Once <F and m are found, (8) gives the 
stress intensity factor at a specific point on the crack front. 
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2.3 Internal Cracks 


The line-spring model described in the previous section may 
easily be extended to internal cracks such as that shown in Fig. 3. 
In this case the basic integral equations remain the same and again 
are given by (1) and (2). However, there are two crack tips (shown 
as A and B in Fig. 3) which create two different stress intensity 
factors for each cross-section taken perpendicular to direction. 
Previous assumptions regarding the stress intensity factors will 
also remain the same; however, they are now defined seperately, 
i. e. , 


= ^/h [<rg At (s) + mg Ab (s)] , 

(39a) 

= '/h"[<Tg Bt (s) + mg Bb (s)] , 

(39b) 


where g At , g Ab , g Bfc , g Bb are the shape functions corresponding to 
crack tips A and B. Unlike the edge crack case, these shape 
functions are now functions of both L(x^)/h and d/h where d is the 
distance from the center of the crack to the neutral plane. Two way 
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parameterization is needed in order to define the shape functions 
which are valid for any kind of internal crack. Since elliptic 
cracks whose major axis are parallel to the free surface are the 
main concern in the present study, parametriza tion can be simplified 
by observing that the distance between the center of any crack 
obtained by taking cross-sections perpendicular to the major axis 
and the center line of the plate is fixed. In other words, d is 
constant for any elliptic crack whose major axis is parallel to the 
free surface. Because of this property, as in the previous case, we 
may again define the shape functions for each crack as a function of 
only one variable, which is L(x.j)/h, as 


ft* 

«At ’ 

(40a) 

__ n 

SAb IcAiS 1 - 1 

v.| 

(40b) 

g Bt = Vifs Y. ^i 32 ^ ^ » 

v.i 

(40c) 

n 


g B b = ^ ^i 31 ** 1 

(40d) 


Keep in mind that coefficients are different for each value of d. 

The next step is the representation of the energy available for 
fracture in terms of both stress intensity factors and product of 
load-load point displacement. As L increases by dL, the energy 
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increment available for fracture may be expressed as 


d(D-V) = 


1-r 


2 


E 


[k| d(L/2) 


+ K| d(L/2)] 


( 41 ) 


giving 


3 

G = — (U-V) 

3l 



+ Kg) 


(42) 


which replaces (11). The rate of energy available for fracture as 
expressed in terms of the load line displacements and forces rem ain s 
the same and is given by (16). If the same procedure as in the 
previous section is followed, the functions cc(ij) may be found as 


L 

a tt = h_1 J ( 4t + 4t )/2 dL » 

°r L 

«bb = h_1 J ( 4b + sib )/2 dL » 

0 L 

a tb = a bt = 11-1 / (g At g Ab + g Bt g Bb )/2 dL * 


(43a) 

(43b) 

(43c) 


After this point, everything, including the resulting integral 
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equations and the solution method, is exactly the same as explained 
in Section 2.2. After solving the equations, the functions G 1 and G 2 
are again substituted to (27) and (28) to find and m, which are 
then substituted to (39) to find the two stress intensity factors at 
the corresponding value of x. 


2.4 Symmetric Internal Crack 

This is a special case of an internal crack where d=0. Due to 
the fact that simple tension cannot create any crack surface 
rotation and simple bending cannot create any crack opening 
displacement Con the neutral plane), from (23) and (24) one can 
write 


«tb = «bt = 0 ’ (W 


which gives 


'*'tb = Y bt = ° ' (2,5) 


Substitution of (45) into (37) would decouple the integral equations 
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(only unknown in (37a) is while the only unknown in (37b) is G 2 ) 
which would in turn enable us to solve them seperately. This 
property reduces both the computer time and memory space needed in 
the numerical solution. 


2.5 Two Opposite Elliptic Edge Cracks 

Two opposite edge cracks shown in Fig. 4 has exactly the same 
properties as symmetric internal cracks except for the coefficients 
of the shape functions. Once the shape functions are established by 
a suitable curve fitting, the same procedure that is applied to 
symmetric internal cracks is applied to find the stress intensity 
factors. 


2.6 Interaction Between Two Identical Cracks 

All the crack geometries that we considered up to now were 
single cracks lying in the region (-a, a) on the x^ axis. Only 
difference among them was their configuration in the thickness 
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direction. Shape functions, energy available for fracture and 
functions cc^ and Y jj are all sensitive to this configuration, 
while the form of the integral equation is independent of it. In 
this section, we will modify (1) and (2) for two identical cracks 
which is shown in Fig. 5. In this case D is defined as (b/a , c/a). 
With the new definition of D, (1) and (2) will take the following 
form. 


C' 


ff(x’) 1 


2<V 




1 1 . 9 

I + j G^fz) dz = , b* < x' < c' , (46) 

*■ z-x ’ z+x 1 J E 


E 


m(x») 


- Ri 


r 1 1 i r 1 1 

+ J- r 3 

L z-x 1 z+x 1 J L ( z-x 1 ' 


( z-x ' ) J ( z+x ' ) 


id 


Rl| [ K 2 (p|z-X'|) + K 2 (p|z+x*|)l 

L z-x' z+x* J 


G 2 (z) dz 


_“o_ 

6E 


, b* < x» < c» , 


(47) 


where 


x' = x./a, b’ = b/a, c' = c/a 


Again we have to convert the limits of integrals to (-1,1) for 
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numerical solution purposes. By introducing the following parameters 
(keep in mind that from Fig. 5, (c’-b'J/B = 1 ) 


x = (2x'-c*-b') / (c’-b* ) , (48) 
x' = [(c'-b'Jx+c'+b’] / 2 = x + [(c'+b*) / 2] , (49) 
dx ' = dx , 

and 

t = (2z-c’-b' ) / (c'-b») , (50) 
z = [ (c’-b* )t+c'+b’] / 2 = t + [(c'+b») / 2] , (51) 
dz = dt , 

and substituting (27) and (28) for ff(x) and m(x) respectively, (46) 
and (47) can be rewritten as follows: 


x 

Y bb (x) / Gf (t ) dt - 


r 2 [— + — ] - h 3 [— + — ] 


^2 


yf 72 


r 4 [ — K 2 (p|y 1 |) +— K 2 (p[y 2 |)1 

h y 2 


G-j ( t ) dt + Y bt (x) G 2 (t) dt 


; 

I 
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-1 < X < 1 


(52) 




6E 




G^t) dt + f tt (x) G 2 (t) dt - 


« 

J 


-I 


1 

2« 

-I 


f [-*- 

Ly 1 y 2 


= , -1 < X < 1 , 

E 


where 


y^i = t - x , 

y 2 =t + x+ b , +c' 


Integral equations (52) and (53) will again be solved 
Chebyshev closed fora integration formula with the 
conditions 


i 

J G^t) dt = 0 , 

J G^t) dt = 0 , 

-l 

and the same backward procedure defined in Section 2,2 
find stress intensity factors at any point on the crack 


] G 2 (t) dt 
(53) 


(54) 

by Gauss- 
additional 


is used to 
defined in 
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the region (b,e). Due to symmetry with respect to the x 2 x^ plane, 
stress intensity factors on the two cracks are also symmetric with 
respect to the same plane. Integral equations are valid for any two 
identical crack (edge, internal, two opposite) problem providing the 
corresponding *f functions defined in Section 2.2 . 


2.7 Interaction of Three Identical Cracks 

As stated earlier integral equations for general case given in 
( 1 ) and ( 2 ) are valid for any number of cracks provided there i 3 
symmetry with respect to the * 2 X 3 plane. This means that any three 
crack system composed of a symmetric crack surrounded, by two 
identical cracks (not needed to be identical with the third one) can 
be solved by using (1) and (2). But in this study, for 
simplification purposes, all of them are chosen identical. 
Configuration and parameters are shown in Fig. 6 . D is defined as 
(0,1),' (b/a , c/a). Because of the nature of the problem, we expect 
that crack surface displacements and rotations be symmetric with 
respect to the x 2 x^ plane. This suggests that two outer cracks 
should have the same rotation and displacement at the equidistant 
points from the origin, while the third one undergoes completely 
different displacement and rotation. For this reason G functions 
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should be defined seperately for inner crack and two outer cracks. 
Let 


G.^ : Crack surface rotation derivative for 
the inner crack, (-a, a), 

G 12 i Crack opening displacement derivative for 
the inner crack, (-a, a), 

G 2 i : Crack surface rotation derivative for 
the outer crack, (b,c), 

G 22 : Crack opening displacement derivative for 
the outer crack, (b,c). 


Also ff(x) and m(x) for the inner and outer cracks should be defined 
seperately as follows, 


ff^x) : ff(x) for inner crack, 
m^ (x) : m(x) for inner crack, 
<r 2 (x) : ff(x) for outer crack, 
m^x) : m(x) for outer crack. 

(1) and (2) can be written for inner crack 


E 



1 

2V 


c 


r 
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1 G 12 (z) dz =■ 

z+x* 


, 0 < x* < 1 , 


(55) 


m^x*) 


- Hi 


6E 


R [ + 1 - R [ + -] 

- ^z-x’ z+x' L (z-x')3 (z+x')^ 


r 1 

1 1 1 

K 2 ( p z-x ' ) + 

K 2 (p z+x* ) 

L z-x' 

z+x* 


- R, 


b' 


\\— + —l - b,[-I— * — — 1 

■ L z-x* z+x ,J '^(z-x'J^ (z+x*)^ 


r 1 1 

+ RJ Kg (pj z-x • l ) + Kg(p| z+x' | ) 

L z-x ' z+x * 


G^ 2 (z) dz 


m„ 


= , 0 < x* < 1 , 

6E 


(56) 


and for outer crack 


ff 2 (x') 1 Cr 1 


2 ir J 

o 


[ -♦ -]g 11 (z) 

Z-X * z+x ' 


dz - 


1 


2ir J, 

b 


1 

z-x* 


z+x' 


1 *0 

J G 12 (z) dz = , b* < x* < c' , 


(57) 
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n^x') 


“ R 1 R 2| 


1 1 
+ 


6E j l “z-x' z+x 

o 


[■ + 1 _ R [ + 1 

L z-x' z+x' J L (z-x')^ (z+x') 3J 


1 1 
+ 


+ Rj, f K- ( p | z-x ' 1 ) + KpCpIz+x* |)1 1 G.p(z) 

L z-x' z+x' J J 


dz 



1 1 


+ ] - R f + 1 

-X' Z+X 1- ' ^(z-x')3 (z+x')^ 


r 1 1 1 

flj K 2 (pjz-x'|) + KgCpjz+x * | )J 

*■ z-x ’ z+x ’ 


G >22 ( z) dz 


m o 

, b' < x* < e» , 

6E 


(58) 


where x'= x/a, b'= b/a, c*=c/a. 
From (36) t we may write 


/[-.-] 

t z-x' z+x' 


r 1 


( z) dz = 


z-x 1 


G^Cz) dz , 


(59) 



1 1 


1 1 


+ - R_[ + - ] + R ij[ K 2 (p|z-x'|) 

-x’ z+x 1 L (z-x') J (z+x')- 5 - 1 L z-x' 
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1 


1 


1 





+ K 2 (p|z+x* I ) 

G^ 2 (z) dz = 


z+x* 

> 



-I 


R 

( Z-X ' ) ^ 


1 

+ Rjj KgCpIz-x' I) 

z-x’ 


G 12 (z) dz 


( 60 ) 


Same expressions as (59) and (60) can also be written for G 21 and 
° 22 ' 


If we define 

x» = x for -1 < x* < 1 

z = t for -1 < z < 1 


and use (49) and (51) for b’ < x* < o' and b' < z < c* 
respectively, and substitute (27), (28) for <*(x) and m(x), (55), 
(56), (57), (5 8) will take the following final form. 


x 

/ G 21 


*r bb (x) J G 2i^ fc ^ dt “ R i 

-I 


1 1 


2 [_ + _] _ R 3 [— 


,y i y 2 


y r 




+ R^r — K 2 (p|y 1 |) + — K 2 (p|y 2 |)] 
*1 y 2 


G 21 (t) dt - R 


'J 

-I 


1 

y 3 


R 2- 
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1 1 

fi 3“i + R ii — Vp|y 3 d 
y 3 y 3 


G n (t) dt + y bt (x) J c^ 2 (t) dt 

-I 


m„ 


— , -1 < x < 1 , 

6 E 


( 61 ) 




G 2 i (t) dt + Ytfc(x) 


tt v 


J 

H 


1 f r 1 

G 22 ( t ) dt [ — 

2V J «■ y, 


1 . 1 f 1 <r 

+ — G 22 (fc) dt — G i2 (t) dt = » -1 < X < 1 , 

y 2 2 *J\ y 3 E 


(62) 


?bb (x > 


G n (t) dt - R 1 


r^JL]_ R 3 r-*-l 

y 4 y 5 3 yj y|-* 


! 


r 1 1 t 

R 4 — * 2 ^ yi* ) +— ^(p y 5 ) I 

y 4 y 5 


G 21 (t) dt - ^ | [Rg — 

y 1 


1 1 

n 3 "V r j» — K 2 (p yi ) 
y i yi 


G n (t) dt + y bt (x ) G 12 (t) dt 


/ 


-I 

x 


m. 


o 

6E 


= — f -1 < x < 1 , 


(63) 


tb 


f * i J i 

(x) G 11 (t) dt + Y t t (x ) I g 12 ^^ dt I" — 

-i 2w h Ly 4 
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1 1 

— 1 G 22 (t) dt 

y 5 J 2V 


1 

— G 12 (t) dt = 


-1 < x < 1 , (64) 


where 


yi 

ys 

y 3 

y^ 

ys 


t-x 

t+x+b'+c* 

t-x - [ (b'+c* ) / 2] 
t+x + [ (b’+c 1 ) / 2] 
t-x + [ (b'+c* ) / 2] 


This problem has four integral equations, (61), (62), (6 3), 
(64), with four unknown functions, G^, G^ 2 , G 2 ^ , G 22 , which will be 
solved by Gauss-Chebyshev closed type integration formula under the 
following single- valuedness conditions. 


f G n (t) dt = 0 

f' 

J G 12 (t) dt * 0 

/ G 21 (t) dt = 0 


f G 22 (t ) dt = 0 

-I 


After solving for the unknown functions, the same backward procedure 
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defined in Section 2.2 will be used for each crack seperately in 
order to find stress intensity factors. 

Provided the corresponding *f functions, equations (61), (62), 
(63), (64) can be applied to any type(edge, internal, collinear) of 
identical cracks. The only reason for considering identical cracks 
to demonstrate interaction of three cracks in this section is 
simplicity. Line-spring model can also be applied to any three crack 
system which has two outer identical cracks with another 
symmetric (with respect to the plane) crack in the middle. In 

this case, y functions and crack length will be defined seperately 
for both(outer ones and inner one) cracks. One should keep in mind 
that crack length is a parameter in the coefficients of (61), (62), 
(63), (64). Attention must be paid to redifine the coefficients in 
such case. 
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3. RESULTS AND DISCUSSION 


There are two main objectives in this study: 

1. Assess the accuracy of line spring model, 

2. Provide extensive numerical results which may be useful 
in application. 

Due to unavailability of the solutions for the cases that we 
are interested in, the first objective can only be reached for the 
central (d/h=0) elliptic internal crack under pure tension. 
Extensive results for central, eccentric internal cracks and two 
opposite edge cracks as well as the interaction among two and three 
cracks are tabulated. In all cases it was assumed that P=0.3 . 
However, the effect of 9 on the stress intensity factors does not 
seem to be significant. 

First step for the application of line spring model is to 
represent the SIFs as a polynomial in L/h. For this purpose, 
extensive information about plane strain (a/L Q ) stress intensity 
factors are needed. These results are obtained from previous 
materials for internal cracks, from [18] for two opposite edge 
cracks and are tabulated in Tables 1,3.5,7,9,11,13,15. Corresponding 
coefficients of shape functions are tabulated in Tables 
2,4,6,8,10,12,14,16. Coefficients of shape functions for the edge 
crack case are directly taken from [ 1 ] and tabulated in Table 17. 
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As noted before, for the application of the line spring model, 
the contour of the part-through crack can be any reasonable curve. 
Elliptic cracks are studied here since it is believed that ellipse 
is the closest contour for the actual shape of the crack which may 
be encountered in practical applications. Thus, crack length for any 
cross section is defined by, 



(x^/a) 2 


L 0 Vl - x 2 


L Q being the total crack length at the midsection (x=0) . Note that 

» 

the limiting values of the SIF are 

K— 0 for a/h— 0 ; K-*>K for a/h — 00 

00 

The first case that is studied is the central internal elliptic 

crack (d/h=0). Extensive results and formulas developed from a 

finite element method for this case are given in [19]. Table 18 and 

Figs. 19,20 show the comparision of the SIFs obtained from this 

study with those generated from the formulas given in [19]. SIFs 

obtained from [19] are represented as K r while the normalizing SIF 
f oo 

is K 0 = <r 0 Y'irL 0 /2 where flr Q = N 22 /h. Both values and the percent 
differences between them are given for various parameters at x=0 and 
x=1/2 . As expected, the SIFs K(L 0 ) for the elliptic crack are 
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consistently smaller than the plane straxn values. The table shows 
that with the exception of relatively small values of a/L 0 at small 
L Q /h (for which the line spring is really not a suitable model) the 
agreement is surprisingly good. From the results, one may conclude 
that line spring can be used with some confidence for any central 
internal elliptic crack which has an a/L 0 ratio larger than 2. 

SIFs for the same geometry, but under pure bending are given in 
Table 20 (Fig. 21). This time, normalizing SIF is defined as 
K Q = m 0 (1TL 0 /2)^ 2 where m Q = 6M^/h . The results given are for the 
tension side. On the compression side the stress intensity factors 
have the same values with a negative sign. Under pure bending, since 
the crack faces on the compression side of the plate would close, 
the. results given in the table cannot be used seperately. The 
results are, of course, useful and valid if the plate is subjected 
to tension, as well as bending, in such a way that the superimposed 
SIF is positive everywhere. 

From Table 20, by excluding the case of a/L Q =0.5 which proved 
to be unreliable, it may be concluded that SIF are independent of 
a/L Q ratio for L Q /h < 0.5 . For L Q /h > 0.5 , SIF begins to increase 
with increasing L Q /h ratio, which is expected. Though we are unable 
to confirm the accuracy of the line spring model in this case, we 
could say that there is no unusual behaviour in the results. Results 
are asymptotically approaching to plane strain values as a/L Q is 


38 



increased except for large L Q /h values. 

Another comparision with the previous finite element results 
[20] is shown in Table 19 (Fig. 23). It should be noted that in the 
results given in Table 19 a/L 0 =1.25 is relatively small for the line 
spring model to be effective. Despite that, the relative error does 
not seem to be very high. 

After symmetric internal cracks, nonsymmetrically located 
internal cracks are studied both under pure tension and pure bending 
conditions. Extensive results which scan almost all possible 
configurations arg given in Tables 21-32 (Figs. 7-18). Normalizing 
SIFs are the same as that defined in the previous case for both 
tension and bending. Again we are unable to verify in which region 
the line spring results are reliable because of the lack of previous 
work done on this kind of problem. Asymptotic convergence to plane 
strain results is fairly good. A close observation of the tension 
results would show that for small a/L 0 ratios SIFs first decrease, 
then increase with increasing L Q /h. This is mainly due to the fact 
that tabulated results are normalized with respect to <r 0 /wL^7i 
Thus, even though the normalized SIF is decreasing, the real SIF may 
increase (which indeed happens in this case) but with a lower rate 
than the rate of increase of L Q . Also, again for small L Q /h and 
a/L Q values SIF on the inside part of the crack front turns out to 
be greater than that on the other part of the crack border which is 
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closer to the free surface. 


Two opposite elliptic edge cracks have been studied only under 
pure tension. Normalizing stress intensity factor is again 

(vL q /2) u c with <F Q = N 22 /h • Decreasing-increasing behaviour of 
the SIFs is again observed here. Same explanation which is given 
earlier is also valid for this case. Plane strain results are very 
close to each other in the region that we observe the slight 
decrease in SIFs. For this reason, a small uncertainity in the 
calculations can easily result with a slight decrease while we are 
expecting a slight increase. Except this behaviour, everything is as 
expected in Table 33 (Fig. 22). 

Results for the interaction of two identical elliptic internal 
cracks under pure bending and pure tension have been tabulated in 
Tables 34 and 35 respectively. It can be seen that there is almost 
no interaction between two identical elliptic edge cracks if the 
distance between them (Fig. 5) is larger than 8 times the half crack 
length, a. Another result which may be observed is that the 
interaction is more effective for the cracks having large L Q /h 
ratios. 

SIFs for the two identical central internal elliptic cracks are 
tabulated in Table 36. Results are given only for tension case 
because it has been observed that there is no appreciable 
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interaction for the crack geometries given in Table 36 under pure 
bending. SIF deviates not more than Ijt from the single crack value 
for b/a=0.1 in bending. It is seen that as the distance between the 
two cracks is increased the single crack solution is easily 
recovered. There is almost no interaction for b/a ^.0.1 for the 
cracks L Q /h < 0.5 . If L Q /h > 0.5 , no interaction region can be 
defined for b/a } 4 . 

Interaction results for any two identical elliptic internal 
cracks are given in Table 37. Seme erratic results are observed for 
small a/L Q ratios at small L Q /h ratios. Clearly, for these crack 
geometries line spring is not a suitable model. 

As the last example, SIFs for the three identical internal 
elliptic cracks under pure tension are calculated and are given in 
Table 38. The only conclusion one may draw from these results is 
that the SIFs on the middle crack are slightly higher than the SIFs 
on the outer crack. But the difference is so small that they may be 
regarded as equal. 
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TABLES 
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Table 1. Stress intensity factors for centrally cracked 
plate subjected to tension N or bending M 
under plane strain conditions. 


L/h 

(K/K 0 ) s 


0.05 

0.1 

1.0060 

0.0250 

0.0500 

0.2 

1.0246 

0.1001 

0.3 

1.0577 

0.1505 

0.4 

1.1094 

0.2023 

0.5 

1.1867 

0.2573 

0.6 

1.3033 

0.3197 

0.7 

1.4884 

0.3986 

0.8 

1.8169 

0.5186 

0.9 

2.585 

0.7776 

0.95 

4.252 

1.1421 


Table 2. The coefficients b^ and c^ for the shape 

- functions g^ and g^ for symmetric internal 
(d/h=0) crack. 


i 

b i 

c i 

1 

0.7070 

0.0169 

2 

0.4325 

-0.4629 

3 

-0.1091 

15.0622 

4 

7.3711 

-143.7384 

5 

-57.7 894 

807.244 9 

6 

271.1551 

-2844 . 852 5 

7 

-744.4204 

6468. 9152 

8 

1183.9529 

-9477.5512 

9 

-1001.4920 

8638.7826 

10 

347.9786 

-4455.2167 

11 


993.24 82 
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Table 3. Stress intensity factors for the plate which has an 
internal crack with d/h=0.05 subjected to uniform 
tension N or bending M under plane strain conditions. 


L/h 

<Wn 

<Wn 

<Wm 

<V k o>M 

0.0001 



0.1001 

0.1000 

0.09 

1.0053 

1.0050 

0.1455 

0.0555 

0.18 

1.0221 

1.0200 

0.1923 

0.0120 

0.27 

1.0530 

1.0452 

0.2409 

-0.0305 

0.36 

1.1021 

1.0823 

0.2 925 

-0.0722 

0.45 

1.1769 

1.1340 

0.3497 

-0.1134 

0.54 

1.2909 

1.2053 

0.4175 

-0.1547 

0.63 

1.4731 

1.3053 

0.5068 

-0.1973 

0.72 

1.7983 

1.4518 

0.6467 

-0.2430 

0.81 

2.5631 

1 .6887 

0.9525 

-0.2923 

0.855 

3.6610 

1.885 8 

1.3814 

-0.3125 


Table 4. The coefficients of b^ and c^ for the shape 
functions g t and g b for internal crack with 
d/h=0.05. 


i 

1 

b Ai 

b Bi 

•H 

O 

c Bi 

1 

0.7071 

0.7071 

0.0708 

0.0707 

2 

0.45 97 

0.4347 

-0.0623 

-0.3701 

3 

0.7671 

-0.0915 

13.1229 

0.5654 

4 

0.1552 

2.6 97 3 

-166.42 80 

-6.6423 

5 

-9.3017 

-14.1195 

1145.8217 

45.7189 

6- 

97.3172 

54. 9653 

-4762.0 914 

-189.9515 

7 

-413.9673 

-135.3432 

12511.5152 

498. 8463 

8 

936.4719 

205.3051 

-20927.0019 

-834.5704 

9 

-1078.2322 

-173-3480 

21613.9362 

862.1672 

10 

504.0555 

62.8847 

-12568.0268 

-501.4354 

11 



3148.4 8T 9 

125.5869 


Table 5. Stress intensity factors for the plate which has an 
internal crack with d/h=0.10 subjected to uniform 
tension N or bending M under plane strain conditions. 


(K a /K 0 ) n (K b /K c ) n (K a /K c ) m (K b /K q ) m 


L/h 


0.0001 

0.08 

0.16 

0.24 

0.32 

0.40 

0.48 

0.56 

0.64 

0.72 

0.76 


1.0050 
1.0212 
1.0513 
1.0998 
1.1743 
1.2887 
1.4722 
1.800 2 
2.5705 
3.6693 


1.0046 
1.0179 
1.0399 
1.0709 
1.1126 
1.1677 
1.241 8 
1.3465 
1.5167 
1.6727 


0.2001 

0.2410 

0.2843 

0.3308 

0.3821 

0.4410 

0.5136 

0.6123 

0.7701 

1.1183 

1.6050 


0.2000 
0.1609 
0.1236 
0.0881 
0.0545 
0.0231 
-0.005 9 
-0.0313 
-0.0504 
-0.0519 
-0.02 80 


Table 6. The coefficients of b A and c ± for the shape 
functions g t and g b for internal crack with 
d/h=0. 10. 


i 

b Ai 

hBi 

c Ai 

c Bi 

1 

0.7071 

0.7072 

0.1415 

0.1414 

2 

'0.5498 

0.5043 

-0. 1734 

-0.3871 

3 

1.5235 

-0.5779 

18.7434 

1.2 936 

4 

-2.2395 

7.64 80 

-266.7713 

-17.0715 

5 

-5.2844 

-52.8793 

2066.46 92 

132.0282 

6' 

226.0267 

257.2074 

-9661.5218 

-617.3023 

7 

-1423.2887 

-799.7410 

28556.2764 

1826.3 1 91 

8 

4348.1446 

1530.8314 

-53734.1216 

-3441.9797 

9 

-6553.5540 

-1634.0240 

62435.9340 

4007.6642 

10 

11 

3959.2116 

749.0673 

-40844.2364 

11511.5912 

-2628.6642 

743.3343 


Table 7. Stress intensity factors for the plate which has an 
internal crack with d/h=0.15 subjected to uniform 
tension N or bending M under plane strain conditions. 


L/h 


<Wn 

< k B /k o>M 

<Wm 

<W 

0.0001 

0.07 

1.0049 

1.0045 

0.3001 

0.3365 

0.3000 

0.2663 

0.14 


1.0208 

1.0172 

0.3763 

0.2352 

0.21 


1.0506 

1.0380 

0.4206 

0.2066 

0.28 


1.0988 

1.0671 

0.4714 

0.1807 

0.35 


1.1729 

1.1057 

0.5319 

0.1580 

0.42 


1.2868 

1.1561 

0.6089 

0.1395 

0.49 


1.46 93 

1.2235 

0.7164 

0.1271 

0.56 


1.7 938 

1.3197 

0.8909 

0.1259 

0.63 


2.5476 

1.4 82 9 

1.2760 

0.1533 

0.665 

3.6065 

1.6426 

1.80 90 

0.2032 


Table 8. The coefficients of b^^ and c ± for the shape 
functions g t and g b for internal crack with 
d/h=0. 15. 


i 

•H 

<< 

JO 

b Bi 

c Ai 

c Bi 

1 

0.7071 

0.7072 

0.2122 

0.2121 

2 

0.7028 

0.6376 

-0.2929 

-0.4042 

3 

2.7653 

-1.2331 

26.3239 

• 2.2494 

4 

-7.2036 

19.00 57 

-427.2558 

-33.6757 

5 

9.1384 

-173.8407 

37 82.95 91 

297.4990 

6 

667.4954 

1108.9410 

-20214.1250 

-1590.1109 

7 

-6105.7233 

-4517.1019 

68285.53 44 

5378.6049 

8 

25260.2847 

11317.3469 

-146859.5866 

-11588.2217 

9 

-50586.0954 

-15802.5485 

195038.2341 

15425.1699 

10 

11 

40325. 8388 

9475.74 80 

-145 833-622 8 
46 980.5243 

-11 566.82 8 8 
3740.3538 


Table 9. Stress Intensity factors for the plate which has an 
internal crack with d/h=0.20 subjected to uniform 
tension N or bending M under plane strain conditions. 


L/h 

<Wn 

<V s o>N 

<Wm 

<W k o>M 

0.0001 



0.4001 

0.4000 

0.06 

1.0048 

1.0044 

0.4319 

0.3718 

0.12 

1.0205 

1.0170 

0.4682 

0.3468 

0.18 

1.0497 

1.0374 

0.5102 

0.3252 

0.24 

1.0968 

1.0660 

0.5601 

0.3070 

0.30 

1.1691 

1.1034 

0.6216 

0.2 930 

0.36 

1.27 99 

1.1535 

0.7020 

0.2844 

0.42 

1.4562 

1.2202 

0.8162 

0.2839 

0.48 

1.7668 

1.3161 

1.0031 

0.2981 

0.54 

2.4756 

1.4 806 

1.4127 

0.3482 

0.57 

3.4498 

1.6422 

1.9674 

0.4165 


Table 10. The coefficients of b^^ and c^ for the shape 
functions g t and g b for internal crack with 
d/h=0.20. 


i 

b Ai 

hBi 

c Ai 

c Bi 

1 

0.7071 

0.7072 

0.2829 

0.2828 

2 

0. 9394 

0.8534 

-0.4105 

-0.4192 

3 

5.0186 

-2.2518 

36.3675 

3.4524 

4 

-19.6345 

47.2610 

-686.9924 

-59.384 8 

5 

76.1489 

-589.3736 

7097.1745 

611.7194 

6 

2376. 8770 

5125.5432 

-44245.1037 

-3814.7743 

7 

-32402.0663 

-2841 3.7 1 81 

174386.6733 

15058.0019 

8 

187563.3073 

96818.06 6 4 

-437594.1661 

-37860.22 90 

9 

-517758.7465 

-183743.9141 

6 7 8087.6 506 

58816.9514 

10 

566112.6482 

14 9736 .5141 

-5 91607.7634 

-51479.6217 

11 



222394.8277 

19433.7456 



Table 11. Stress intensity factors for the plate which has an 
internal crack with d/h=0.25 subjected to uniform 
tension N or bending M under plane strain conditions. 


L/h 


<Wn 

<Wm 

<Wm 

0.0001 



0.5001 

0.5000 

0.05 

1.0046 

1.0042 

0.5273 

0.4771 

0.10 

1.0197 

1.0165 

0.5598 

0.45 83 

0.15 

1.0476 

1.0364 

0.5991 

0.4434 

0.20 

1.0925 

1.0644 

0.6473 

0 . 4327 

0.25 

1.1610 

1.1016 

0.7084 

0.4270 

0.30 

1.2652 

1.1505 

0.7900 

0.4277 

0.35 

1.4295 

1.2160 

0.9074 

0.4380 

0.40 

1.7151 

1.3099 

1.0996 

0.464 9 

0.45 

2.3529 

1.46 84 

1.5149 

0.5310 

0.475 

3.2077 

1.6192 

2.0658 

0.6085 


Table 12. The coefficients of b A and ^ for the shape 
functions g fc and g b for internal crack with 
d/h=0.25. 


i 

•H 

<«: 

XI 

b Bi 

c Ai 

c Bi 

1 

0.7071 

0.7071 

0.3536 

0.3536 

2 

1.3041 

1 . 1 91 8 

-0.52 46 

-0.4284 

3 

9.6510 

-4.17 93 

50.2589 

4.8592 

4 

-57.8163 

129.7358 

-1136.1027 

-97.65 5 8 

5 

450.7761 

-2325.3551 

14085.0245 

1206.5210 

8 

10351.8621 

29069.6053 

-105370.0173 

-9029. 9162 

7 

-227973.9193 

-231101.0011 

4 98386.147 5 

42788.3389 

8 

1 960136.6080 

1128028.3541 

-1500844.1305 

-12 9153.3052 

9 

-7899243.5583 

-3063223.0225 

27 91092.0 4 97 

240895.9775 

10 

11 

12538854.7550 

3568834.50 97 

-2922538.277 9 
13185 91.9416 

-253 1 69.5873 
114779.2965 
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Table 13. Stress intensity factors for the plate which has an 
internal crack with d/h=0.30 subjected to uniform 
tension N or bending M under plane strain conditions. 


W h 



<Wm 

<Wm 

0.001 



0.6001 

0.6000 

0.04 

1.0043 

1.0040 

0.6226 

0.5824 

0.08 

1.0183 

1.0155 

0.6510 

0. 5686 

0.12 

1.0442 

1.0343 

0.6867 

0.5607 

0.16 

1.0855 

1.0608 

0.7320 

0.5568 

0.20 

1.1481 

1.0959 

0.7909 

0.5584 

0.24 

1.2425 

1.1419 

0.8708 

0.5669 

0.28 

1.3898 

1.2031 

0.9866 

0.5855 

0.32 

1.6423 

1.2895 

1.1755 

0.6209 

0.36 

2.1937 

1.4309 

1.5773 

0.6939 

0.38 

2.9157 

1.55 91 

2.0992 

0.7709 


Table 14. The coefficients of b^ and c^ for the shape 
functions g t and g b for internal crack with 
d/hr0.30. 


b Ai 


1j 0.7071 

1.9027 
20.8636 
-197.8895 
2675.1513 
- 71601.3880 
-2639273.5314 
35994016.1511 
227 958779.21 55 

10 567162515.8604 

11 


hi 


0.7071 
1. 7*180 
-8. 9087 
440.7103 
-123 90.96 6 8 
241718.2677 
-2989605.6659 
226 1 8755. 1306 
-94950523 . 9757 
170542044.2578 


c Ai 


0.4244 
-0.6553 
72.3179 
-2037.1553 
3 1 56 9.7 5 0 5 
-295202.5115 
17 4 5346.790 8 
-6570139-7052 
1527381 8.2302 
■19993047.6604 
1127 6 922.90 53 


c Bi 


0.4243 
-0.4321 
6.6403 
-160.0429 
2468.7747 
-23 0 90.84 86 
136788.4159 
-516213.4074 
120392 8.2026 
-1582260.061 4 
897 2 40.7 9 83 


Table 15. "Stress intensity factors for the plate with 
collinear edge cracks subjected to uniform 
tension N under plane strain conditions. 


L/h 

«/Vn 

0.0001 

1.1221 

0.1 

1.1231 

0.2 

1.1254 

0.3 

1.12 92 

0.4 

1.1370 

0.5 

1.1546 

0.6 

1.2117 

0.7 

1.3254 

0.8 

1.5393 

0.9 

2.0836 


Table 16. The coefficients of b^. for toe shape function 
for collinear edge crack. 


i 

b i 

1 

0.7 934 

2 

0.0775 

3 

-0.7542 

4 

7.5825 

5 

-12. 1712 

6 

-186.5011 

7 

1236.2 85 8 

8 

-3043.6190 

9 

3350.3456 

10 

-1374.8426 
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Table 17. The coefficients of and c ± for the shape 
functions g t and gjj for edge crack. 


i 

b i 

c i 

1 

1.1216 

1.1202 

2 

6.5200 

-1.8872 

3 

-12.3877 

18.0143 


89.0554 

-87.3851 

5 

-188.6080 

241.9124 

6 

207.3870 

-319.9402 

7 

-32.0524 

168.0105 
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Table 18. Comparison of the stress intensity factors K(x) 

calculated in this study at x=0 and x= 1/2 (xsx^/a) 
for an internal planar elliptic crack in a plate under 
uniform tension N with the corresponding values K r (x) 
given in Ref. [ 19 ] . ?D=100(K r -K)/K r . 


L 0 /h 

a/L 0 

X 

K(x)/K 0 

K r (x)/K 0 

$D 

0.1 

0.5 

0 

0 . 916 

0.637 

43.7 

0.1 

0.5 

1/2 

0.86 8 

0.637 

36.2 

0.1 

1.0 

0 

0.955 

0.827 

15.5 

0.1 

1.0 

1/2 

0.896 

0.785 

14.2 

0.1 

2.0 

0 

0.976 

0.935 

4.3 

0.1 

2.0 

1/2 

0.911 

0.875 

4.1 

0.1 

3.0 

0 

0.983 

0.967 

1.7 

0.1 

3.0 

1/2 

0.916 

0.902 

1.6 

0.1 

4.0 

0 

0.987 

0.980 

0.6 

0.1 

4.0 

1/2 

0.919 

0.914 

0.5 

0. 1 

10.0 

0 

0.993 

0.999 

- 0.6 

0.1 

10.0 

1/2 

0.923 

0 . 930 

- 0.7 

0.2 

0.5 

0 

0 . 862 

0.638 

35. 1 

0.2 

0.5 

1/2 

0.827 

0.638 

2 9.6 

0.2 

1.0 

0 

0.931 

0 . 830 

12.2 

0.2 

1.0 

1/2 

0.880 

0.788 

11.6 

0.2 

2.0 

0 

0.971 

0.942 

3 . 1 

. 0.2 

2.0 

1/2 

0.90 8 

0.881 

3.1 

0.2 

3.0 

0 

0.986 

0 . 976 

1.0 

0.2 

3.0 

1/2 

0.91 8 

0.910 

0.9 

0.2 

4.0 

0 

0.993 

0.991 

0.2 

0.2 

4.0 

1/2 

0.923 

0.923 

0.0 

0.2 

10.0 

0 

1.007 

1.013 

- 0.6 

0.2 

10.0 

1/2 

0.933 

0.942 

- 1.0 
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Table 18 - Cont 


V h 

a / L 0 

X 

K ( x )/ K 0 

K r ( x )/ K 0 

?D 

0.3 

0.5 

0 

0.824 

0.641 

28.6 

0.3 

0.5 

1/2 

0.796 

0.640 

24.4 

0.3 

1.0 

0 

0.920 

0.837 

9.9 

0.3 

1.0 

1/2 

0.871 

0.794 

9.8 

0.3 

2.0 

0 

0.979 

0.957 

2.3 

0.3 

2.0 

1/2 

0.914 

0.893 

2.3 

0.3 

3.0 

0 

1.001 

0.996 

0.5 

0.3 

3.0 

1/2 

0.930 

0.927 

0.3 

0.3 

4.0 

0 

1.012 

1 . 014 . 

- 0.2 

0.3 

4.0 

1/2 

0.937 

0.942 

- 0.5 

0.3 

10.0 

0 

1.034 

1.041 

- 0.7 

0.3 

10.0 

1/2 

0.952 

0.966 

- 1.5 

0.4 

0.5 

0 

0.798 

0.645 

23.6 

0.4 

0.5 

1/2 

0.775 

0.644 

20.3 

0.4 

1.0 

0 

0.920 

0.851 

8. 1 

0.4 

1.0 

1/2 

0.871 

0.804 

8.3 

0.4 

2.0 

0 

1.000 

0.984 

1.6 

0.4 

2.0 

1/2 

0.92 9 

0.915 

1.6 

' 0.4 

3.0 

0 

1.030 

1.031 

-0. 1 

0.4 

3.0 

1/2 

0.950 

0.955 

- 0.4 

0.4 

4.0 

0 

1.047 

1.054 

- 0.7 

0.4 

4.0 

1/2 

0.961 

0.974 

- 1.3 

0.4 

10.0 

0 

1. 078 

1.0 91 

- 1.2 

0.4 

10.0 

1/2 

0.982 

1.005 

- 2.3 
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Table 18 - Cont 


L 0 /h 

a / L 0 

X 

K ( x )/ K 0 

K r ( x )/ K 0 

*D 

0.5 

0.5 

0 

0.783 

0.654 

19.8 

0.5 

0.5 

1/2 

0.761 

0.650 

17.2 

0.5 

1.0 

0 

0 . 932 

0.874 

6.7 

0.5 

1.0 

1/2 

0.880 

0.821 

7.3 

0.5 

2.0 

0 

1.036 

1.030 

0.6 

0.5 

2.0 

1/2 

0.956 

0.949 

0.7 

0.5 

3.0 

0 

1.078 

1.090 

-1.1 

0.5 

3.0 

1/2 

0.984 

0.998 

- 1.4 

0.5 

4.0 

0 

1.101 

1.121 

- 1.8 

0.5 

4.0 

1/2 

0.998 

1.023 

- 2.4 

0.5 

10.0 

0 

1.145 

1.172 

- 2.4 

0.5 

10.0 

1/2 

1.025 

1.063 

- 3.5 

0.6 

0.5 

0 

0.779 

0.667 

16 . 9 

0.6 

0.5 

1/2 

0.756 

0.658 

14.9 

0.6 

1.0 

0 

0 . 960 

0 . 911 

5.4 

0.6 

1.0 

1/2 

0.901 

0.846 

6.4 

0.6 

2.0 

0 

1.095 

1.103 

- 0.8 

0.6 

2.0 

1/2 

0.997 

0.999 

- 0.2 

0.6 

3.0 

0 

1.151 

1. 183 

- 2.7 

0.6 

3.0 

1/2 

1.033 

1.061 

- 2.6 

0.6 

4.0 

0 

1.183 

1.225 

- 3-5 

0.6 

4.0 

1/2 

1.052 

1.092 

- 3.6 

0.6 

10.0 

0 

1.245 

1.298 

- 4.0 

0.6 

10.0 

1/2 

1.088 

1.143 

- 4.7 
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Table 18 - Cont 


L 0 /h 

a / L 0 

X 

K ( x )/ K 0 

K r ( x )/ K 0 

?D 

0.7 

0.5 

0 

0.788 

0.687 

14.7 

0.7 

0.5 

1/2 

0.760 

0.671 

13.3 

0.7 

1.0 

0 

1.009 

0 . 968 

4.3 

0.7 

1.0 

1/2 

0.935 

0.882 

6.0 

0.7 

2.0 

0 

1.187 

1.213 

- 2.2 

0.7 

2.0 

1/2 

1.058 

1.067 

- 0.8 

0.7 

3.0 

0 

1.266 

1.322 

- 4.2 

0.7 

3.0 

1/2 

1.106 

1.144 

- 3.3 

0.7 

4.0 

0 

1.310 

1.381 

- 5. 1 

0.7 

4.0 

1/2 

1.132 

1.183 

- 4.4 

0.7 

10.0 

0 

1.403 

1.4 83 

- 5.4 

0.7 

10.0 

1/2 

1.179 

1.242 

- 5.1 

0 . 8 

0.5 

0 

0 . 818 

0.717 

14. 1 

0.8 

0.5 

1/2 

0.777 

0.686 

12.9 

0.8 

1.0 

0 

1.096 

1.051 

4.3 

0.8 

1.0 

1/2 

0.991 

0.930 

6.6 

0 . 8 

2.0 

0 

1.341 

1.372 

- 2.3 

0.8 

2.0 

1/2 

1.152 

1.152 

0.0 

0.8 

3.0 

o 

1.457 

1.521 

- 4.3 

0.8 

3.0 

1/2 

1.218 

1.245 

- 2.2 

0.8 

4.0 

0 

1.525 

1.603 

— 4.9 

0.8 

4.0 

1/2 

1.253 

1.290 

- 2.9 

0 . 8 

10.0 

0 

1.674 

1.747 

- 4.2 

0.8 

10.0 

1/2 

1.318 

1.349 

- 2.3 
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Table 1 8- Cont 


L c /h 

a/L 

o 

X 

K(x)/K q 

K r (x)/K 0 

ZD 

0.9 

0.5 

0 

0.905 

0.75 9 

19.3 

0.9 

0.5 

1/2 

0.813 

0.709 

14.7 

0.9 

1.0 

0 

1.284 

1.167 

10.0 

0.9 

1.0 

1/2 

1.086 

0.987 

10.1 

0.9 

2.0 

0 

1.661 

1.594 

4.2 

0.9 

2.0 

1/2 

1.310 

1.246 

5.1 

0.9 

3.0 

0 

1.858 

1.799 

3.3 

0.9 

3.0 

1/2 

1.405 

1.347 

4.2 

0.9 

4.0 

0 

1.981 

1.912 

3.6 

0.9 

4.0 

1/2 

1.456 

1.392 

4.6 


Table 19. Comparison of the stress intensity factors K(L Q ) 
calculated in this study at the midsection of a 
symmetrically located internal (d/h=0) planar crack 
in a plate under uniform tension N with the 
corresponding results K r given in [20]. L /h=.75, 
a/L 0 = 1.25, x=x 1 /a=cos8, *D = 1 00(K-K r )/K r . 


e 

X 

K r /K 0 

K/K 0 

ZD 

90° 

0 

0.985 

1.120 

13.7 

80° 

0. 174 

0. 971 

1.103 

13.6 

70° 

0.342 

0.944 

1.052 

11.4 

pa 

0.500 

0.898 

0.S73 

8.4 

mSm 

0.707 

0.810 

0.832 

2.7 

40° 

0.766 

0.770 

0.742 

-3.6 
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Table 20. Normalized stress intensity factors calculated 
at the midsection of a symmetrically located 
(d/h=0) internal planar elliptic crack in a 
plate under uniform bending M. 


L 0 /h ‘ 

l/L 0 -» 0.5 

2.0 

3.0 

4.0 

5.0 

10.0 

Plane 

Strain 

0.1 

.050 

.050 

.050 

.050 

.050 

.050 

.050 

0.2 

.099 

.099 

.099 

.099 

.099 

.099 

.100 

0.3 

.148 

.149 

.149 

.149 

.149 

.149 

.151 

0.4 

.195 

.199 

.199 

.199 

.200 

.200 

.202 

0.5 

.240 

.250 

.251 

.252 

.253 

.254 

.257 

0.6 

.285 

.305 

.308 

.310 

.311 

.313 

.320 

0.7 

.332 

.369 

.375 

.379 

.382 

.387 

.399 

0.8 

.387 

.457 

.470 

.478 

.4 83 

.4 96 

.519 

0.9 

.476 

.619 

.650 

.670 

.683 

.717 

.778 
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Table 21. Normalized stress intensity factors calculated at the 
midsection of an internal planar elliptic crack with 
d/h=.05 in a plate under uniform tension N. 


V*o 


a 

L 0 /h 

/L o ^0.5 

1.0 

2.0 

3.0 

4.0 

7.0 

10.0 

Plane 

Strain 

.05 

0.954 

0.975 

0.987 

0.990 

0.992 

0.995 

0.996 

1.002 

.15 

0.887 

0.944 

0. 976 

0.988 

0.994 

1.001 

1.004 

1.015 

.25 

0.838 

0.926 

0.979 

0.998 

1.008 

1.022 

1.028 

1.045 

.35 

0.804 

0.921 

0.996 

1.025 

1.040 

1.061 

1.070 

1.096 

.45- 

0.782 

0.931 

1.032 

1.073 

1.095 

1.126 

1.139 

1.177 

.55 

0.775 

0.959 

1.096 

1.153 

1.185 

1.231 

1.251 

1.307 

.65 

0.786 

1.018 

1.204 

1.288 

1.336 

1.407 

1.439 

1.529 

.75 

0. 828 

1.134 

1.407 

1.539 

1.619 

1.742 

1.799 

1.970 

.85 

0.994 

1.479 

1.987 

2.268 

2.452 

2.761 

2.921 

3.481 





v * 0 

a/L n ^0.5 

1.0 

2.0 

V h 




.05 

0.955 

0.976 

0.987 

.15 

0.891 

0.946 

0. 977 

.25 

0.847 

0.92 9 

0.978 

.35 - 

0. 817 

0. 923 

0.991 

.45 

0.798 

0.92 8 

1.017 

.55 

0.789 

0.944 

1.056 

.65 

0.788 

0.972 

1.115 

.75 

0.797 

1.017 

1.201 

.85 

0.822 

1.088 

1.334 


3.0 

4.0 

7.0 

10.0 

Plane 

Strain 

0.991 

0.993 

0.995 

0.996 

1.002 

0.988 

0.993 

1.001 

1.003 

1.014 

0.996 

1.006 

1.018 

1.023 

1.039 

1.017 

1.031 

1.049 

1.056 

1.078 

1.051 

1.070 

1.095 

1.105 

1.134 

1.102 

1.127 

1.161 

1. 176 

1.215 

1.176 

1.210 

1.257 

1.277 

1.333 

1.284 

1.332 

1.401 

1.432 

1.517 

1.455 

1.529 

1.643 

1.697 

1.85 9 
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Table 22. Normalisred stress intensity factors calculated at the 
midsection of an internal planar elliptic crack with 
d/h=0.10 in a plate under uniform tension N. 


k a /k 


o 


L 0 /h 

a/ 

’L o -0.5 

1.0 

2.0 

3.0 

4.0 

7.0 

IQ. 0 

Plane 

Strain 

.05 


0.950 

0.974 

0.986 

0.990 

0.992 

0.995 

0.996 

1.002 

.15 


0. 878 

0.940 

0. 975 

0.988 

0.994 

1.003 

1.006 

1.018 

.25 


0.825 

0.921 

0.980 

1.002 

1.013 

1.029 

1.036 

1.056 

.35 


0.789 

0. 919 

1.004 

1.037 

1.055 

1.081 

1.092 

1.124 

.45 


0.770 

0.937 

1.056 

1.105 

1.132 

1.172 

1.189 

1.240 

.55 


0.770 

0.986 

1.152 

1.226 

1.268 

1.331 

1.360 

1.444 

.65 


0.805 

1.093 

1.342 

1.462 

1.534 

1.646 

1.699 

1.860 

.75 


0. 961 

1.426 

1.900 

2.159 

2.329 

2.615 

2.764 

3.291 





V K o 






T # 

a/L —0.5 

i ^ 

1.0 

2.0 

3.0 

4.0 

7.0 

10.0 

Plane 


L 0 /h 

Strain 

.05 

0.952 

0.974 

0.986 

0.990 

0.993 

0.995 

0.996 

1.002 

.15 

0.886 

0.943 

0. 976 

0.988 

0.993 

1.001 

1.005 

1.016 

.25 

0.843 

0.927 

0.979 

0.998 

1.007 

1.021 

1.026 

1.043 

.35 

0.815 

0.923 

0.993 

1.020 

1.034 

1.054 

1.062 

1.085 

.45 

0.798 

0.92 8 

1.019 

1.055 

1.075 

1.102 

1.113 

1.145 

.55 - 

0.789 

0.944 

1.059 

1.107 

1.133 

1. 171 

1.187 

1.231 

.65 

0.788 

0.970 

1.116 

1.181 

1.218 

1.272 

1.296 

1.363 

.75 

0.795 

1.010 

1.203 

1.297 

1.355 

1.445 

1.488 

1.625 
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Table 23* Normalized stress intensity factors calculated at the 
midsection of an internal planar elliptic crack with 
d/h=0.15 in a plate under uni form tension N. 


k a /k 


o 


a/L 0 -^°.5 1.0 2.0 3.0 4. 0 7.0 10.0 Plane 

L„/h I Straii 


L 0 /h 

Strain 

.05 

0.944 

0.971 

0.985 

0.989 

0.992 

0.995 

0.996 

1.002 

.15 

0.864 

0. 934 

0. 974 

0.988 

0.996 

1.006 

1.010 

1.024 

.25 

0.80 8 

0.916 

0.984 

1.009 

1.022 

1.041 

1.049 

1.076 

.35 

0.773 

0.921 

1.020 

1.060 

1.083 

1.115 

1.129 

1.173 

.45 

0.761 

0.956 

1.101 

1.164 

1.200 

1.254 

1.279 

1.354 

.55 

0.784 

1.051 

1.273 

1.378 

1.441 

1.540 

1.588 

1.733 

.65 

0.926 

1.362 

1.789 

2.021 

2.172 

2.427 

2. 560 

3.033 


k b /k 


o 


a/L -*0.5 1.0 

V h I 


.05 

0.964 

0.972 

.15 

0.&6 

0. 939 

.25 

0.833 

0.923 

.35 

0. 806 

0. 920 

.45 

0.791 

0.92 9 

.55 

0.784 

0.947 

.65 

0.783 

0.976 


2.0 3.0 4.0 


0.985 

0.990 

0.992 

0. 97 5 

0.988 

0.994 

0.979 

1.000 

1.012 

0.997 

1.027 

1.044 

1.028 

1.070 

1.093 

1.076 

1.134 

1.167 

1.152 

1.239 

1.293 


7.0 10.0 Plane 

Strain 


0.995 

0.996 

1.002 

1.003 

1.007 

1.020 

1.027 

1.033 

1.054 

1.067 

1.077 

1.106 

1.126 

1.141 

1.182 

1.216 

1.238 

1.303 

1.379 

1.422 

1.561 
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Table 24. Normalized stress intensity factors calculated at the 
midsection of an internal planar elliptic crack with 
d/h=0.20 in a plate under uniform tension N. 


V K o 


a/ 

L 0 /h 

V-0.5 

1.0 

o 

. 

CM 

3.0 

4.0 

7.0 

10.0 

Plane 

Strain 

.05 

0.964 

0.967 

0.983 

0.989 

0.991 

0.995 

0.996 

1.003 

.15 

0.848 

0.928 

0.974 

0.990 

0.999 

1.011 

1.016 

1.033 

.25 

0.790 

0.914 

0.992 

1.022 

1.039 

1.062 

1.073 

1.107 

.35 

0.761 

0.933 

1.054 

1.104 

1.133 

1.176 

1.196 

1.258 

.45 

0.770 

1.019 

1.202 

1.289 

1.342 

1.424 

1.463 

1.587 

.55 

0.8S6 

1.2 93 

1.666 

1.863 

1.990 

2.204 

2.316 

2.710 





k b /k ( 

3 





a/L_-*.0. 5 

1.0 

2.0 

3.0 

4.0 

7.0 

10.0 

Plane 

L 0 /h 

V 







Strain 

.05 

0.939 

0.968 

0.984 

0.989 

0.992 

0.995 

0.996 

1.003 

.15 

0.863 

0. 933 

0.974 

0.988 

0.996 

1.006 

1.011 

1.026 

.25 

0.81 8 

0.91 8 

0.982 

1.006 

1.020 

1.038 

1.046 

1.072 

.35 

0.793 

0. 918 

1.007 

1.043 

1.063 

1.092 

1.105 

1.144 

.45 

0.779 

0.931 

1.049 

1.102 

1.132 

1.178 

1.199 

1.263 

.55 

0.770 

0.955 

1.121 

1.204 

1.256 

1.340 

1.382 

1.524 
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Table 25. Normalized stress intensity factors calculated at the 
midsection of an internal planar elliptic crack with 
d/h=0.25 in a plate under uniform tension N. 


k a /k 


o 


a/L 0 — 0.5 1.0 2.0 3-0 4.0 7.0 10.0 Plane 


L 0 /h 

Strain 

.05 

0.927 

0.962 

0.981 

0.988 

0.991 

0.995 

0.997 

1.005 

.15 

0.831 

0.922 

0.976 

0.995 

1.006 

1.020 

1.026 

1.04 8 

.25 

0.776 

0.919 

1.013 

1.051 

1.071 

1.101 

1.115 

1.161 

.35 

0.768 

0.977 

1.133 

1.201 

1.241 

1.303 

1.333 

1.42 9 

.45 

0.876 

1.228 

1.540 

1.698 

1.798 

1.964 

2.051 

2.353 


V K 


o 


a/L -*.0.5 

1.0 

2.0 

V h 




.05 

0.930 

0. 964 

0.982 

.15 

0.847 

0.926 

0.973 

.25 

0.801 

0.913 

0.989 

.35 

0.775 

0. 919 

1.026 

.45 

0.758 

0.938 

1.095 


3.0 4.0 7.0 10.0 Plane 

Strain 


0.988 

0.991 

0.995 

0.997 

1.004 

0.991 

1.000 

1.012 

1.017 

1.036 

1.018 

1.034 

1.057 

1.068 

1.102 

1.072 

1.098 

1.138 

1.157 

1.216 

1.172 

1.220 

1.297 

1.336 

1.468 
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Table 26. Normalized stress intensity factors calculated at the 
midsection of an internal planar elliptic crack with 
d/h=0.30 in a plate under uniform tension N. 


k a /k 


o 


a/L o -*0.5 1.0 2.0 3.0 4.0 7.0 10.0 Plane 


L 0 /h 

Strain 

.05 

0.917 

0.958 

0.980 

0.987 

0.991 

0.996 

0.998 

1.007 

.15 

0.816 

0.921 

0.985 

1.008 

1.021 

1.038 

1.046 

1.073 

.25 

0.780 

0.951 

1.070 

1.119 

1.146 

1.187 

1.207 

1.273 

.35 

0. 869 

1.170 

1.420 

1.537 

1.610 

1.727 

1.787 

1.996 


V K o 


a/L rt — 0.5 


1.0 


2.0 


3.0 


4.0 


7.0 10.0 Plane 


L 0 /h 

Strain 

.05 

0.919 

0.959 

0.980 

0.987 

0.991 

0.996 

0.997 

1.006 

.15 

0.829 

0.920 

0. 976 

0.997 

1.008 

1.023 

1.030 

1.053 

.25 

0.783 

0.914 

1.005 

1.042 

1.063 

1.093 

1.108 

1.156 

.35 

0.753 

0.927 

1.071 

1.138 

1. 178 

1.243 

1.276 

1.387 
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Table 27. Normalized stress intensity factors calculated at the 
midsection of an internal planar elliptic crack with 
d/h=0.05 in a plate under uniform bending M. 


k a /k 


o 


a/L„— 1.0 

O 

• 

CM 

3.0 

L 0 /h 




.05 

0.118 

0.11 9 

0.120 

.15 

0.170 

0.173 

0.174 

.25 

0.217 

0.222 

0.224 

.35 

0.266 

0.274 

0.27 8 

.45 

0.316 

0.329 

0.335 

.55 

0.371 

0.391 

0.400 

.65 

0.436 

0.469 

0.484 

.75 

0.528 

0.586 

0.614 

.85 

0.753 

0.889 

0.963 


4.0 7.0 10.0 Plane 

Strain 


0.120 

0.120 

0.120 

0.121 

0.175 

0. 176 

0.176 

0.177 

0.225 

0.227 

0.228 

0.230 

0.279 

0.282 

0.283 

0.287 

0.338 

0.342 

0.344 

0.350 

0.405 

0.412 

0.416 

0.426 

0.4 93 

0.506 

0.513 

0.532 

0.631 

0.659 

0.673 

0.716 

1.013 

1.099 

1.144 

1.313 


V K o 


a/L_— ■ 1 . 0 

2.0 

3.0 

L 0 /h 




.05 

0.073 

0.074 

0.074 

.15 

0.020 

0.023 

0.024 

.25 

-0.032 

-0.027 

-0.025 

.35 

-0.082 

-0.075 

-0.073 

.45 

-0.130 

-0.123 

-0.120 

.55 

-0.178 

-0. 170 

-0.167 

.65 

-0.226 

-0.21 9 

-0.215 

.75 

-0.278 

-0.272 

-0.269 

.85 

-0.339 

-0.338 

-0.334 


4.0 

7.0 

10.0 

Plane 

Strain 

0.074 

0.074 

0.075 

0.075 

0.025 

0.025 

0.026 

0.026 

-0.024 

-0.023 

-0.022 

-0.021 

-0.072 

-0.070 

-0.069 

-0.068 

-0.118 

-0.116 

-0.115 

-0.113 

-0.165 

-0.162 

-0.161 

-0.159 

-0.213 

-0.210 

-0.209 

-0.207 

-0.267 

-0.264 

-0.262 

-0.259 

-0.331 

-0.325 

-0.322 

-0.311 
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Table 28. Normalized stress intensity factors calculated at the 
midsection of an internal planar elliptic crack with 
d/h=0.10 in a plate under uniform bending M. 


V K o 


a/L — 1.0 

O 

• 

(M 

3.0 

4.0 

7.0 

10.0 

Plane 

L 0 /h 







Strain 

.05 

0.216 

0.218 

0.219 

0.219 

0.220 

0.220 

0.221 

.15 

0.263 

0.270 

0.273 

0.274 

0.276 

0.276 

0.279 

.25 

0.308 

0.321 

0.325 

0.328 

0.331 

0.332 

0.337 

.35 

0.357 

0.376 

0.383 

0.387 

0.393 

0.395 

0.403 

.45 

0.412 

0.440 

0.451 

0.457 

0.467 

0.471 

0.484 

.55 

0.477 

0.520 

0.538 

0.549 

0.566 

0.574 

0.598 

.65 

0.571 

0.642 

0.676 

0.697 

0.730 

0.746 

0.797 

.75 

0. 804 

0.959 

1.045 

1.101 

1.199 

1.251 

1.441 

- 



V K o 





a/L — 1 . 0 

2.0 

3.0 

4.0 

7.0 

10.0 

Plane 

V h 







Strain 

.05 

0.170 

0.172 

0.173 

0.173 

0.174 

0.174 

0.175 

.15 

0.114 

0.120 

0.123 

0.124 

0/125 

0.126 

0.128 

.25 

0.061 

0.071 

0.075 

0.077 

0.080 

0.081 

0.084 

.35 

0.010 

0.024 

0.030 

0.033 

0.037 

0.038 

0.043 

.45 

-0.038 

-0.020 

-0.013 

-0.009 

-0.004 

-0.001 

0.005 

.55- 

-0.084 

-0. 062 

-0.053 

-0.048 

-0.040 

-0.037 

-0.028 

.65 

-0.128 

-0.102 

-0.089 

-0.082 

-0.071 

-0.066 

-0.052 

.75 

-0.171 

-0.137 

-0.117 

-0.105 

-0.083 

-0.073 

-0.037 
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Table 29. Normalized stress intensity factors calculated at the 
midsection of an internal planar elliptic crack with 
d/h=0.15 in a plate under uniform bending M. 


K A /K o 


a 

L 0 /h 

'L 0 — 1.0 

2.0 

3.0 

4.0 

7.0 

10.0 

Plane 

Strain 

.05 

0.313 

0.317 

0.319 

0.319 

0.320 

0.321 

0.323 

.15 

0.355 

0.367 

0.371 

0.373 

0.376 

0.378 

0.382 

.25 

0.399 

0.420 

0.428 

0.432 

0.438 

0.440 

0.449 

.35 

0.450 

0.482 

0.495 

0.502 

0.512 

0.517 

0.532 

.45 

0.514 

0.563 

0.564 

0.596 

0.614 

0.623 

0.650 

.55 

0.610 

0.689 

0.727 

0.750 

0.786 

0.804 

0.860 

.65 

0.846 

1.015 

1.106 

1.167 

1.270 

1.325 

1.524 


a/ 

L 0 /h 

'L 0 — 1.0 

2.0 

3.0 

4.0 

7.0 

10.0 

Plane 

Strain 

.05 

0.266 

0.270 

0.272 

0.273 

0.273 

0.274 

0.275 

.15 

0.207 

0.218 

0.221 

0.223 

0.226 

0.227 

0.231 

.25 

0.152 

0.169 

0.176 

0.179 

0.184 

0.1 86 

0.191 

.35 

0.101 

0.125 

0.134 

0.139 

0.146 

0.14 9 

0.158 

.45 

0.055 

0.085 

0.098 

0.105 

0.116 

0.120 

0.133 

.55 

0.011 

0.051 

0.069 

0.080 

0.096 

0.104 

0.125 

.65 

-0.028 

0.028 

0.058 

0.077 

0.108 

0.123 

0.176 
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Table 30. Normalized stress intensity factors calculated at the 
midsection of an internal planar elliptic crack with 
d/h=0.20 in a plate under uniform bending M. 


V K o 


a/L n —*-1.0 

2.0 

3.0 

L 0 /h 




.05 

0.410 

0.417 

0.419 

.15 

0.445 

0. 464 

0.471 

.25 

0.490 

0.522 

0.534 

.35 

0.549 

0.600 

0.621 

.45 

0.642 

0.725 

0.764 

.55 

0. 875 

1.050 

1.142 


4.0 7.0 10.0 Plane 

Strain 


0.420 

0.421 

0.422 

0.425 

0.474 

0. 479 

0.4 81 

0.488 

0.541 

0.551 

0.555 

0.569 

0.633 

0.651 

0.660 

0.687 

0.787 

0.823 

0.841 

0.897 

1.202 

1.303 

1.357 

1.54 9 


k b /k 


o 


a/L — *.1 .0 

2.0 

3.0 

L 0 /h 




.05 

0.362 

0.368 

0.370 

.15 

0.298 

0.315 

0.321 

.25 

0.242 

0.268 

0.278 

.35 

0.193 

0.229 

0.244 

.45 

0.148 

0.197 

0.219 

.55 

0.109 

0.181 

0.218 


4.0 7.0 10.0 Plane 

Strain 


0.372 

0.373 

0.373 

0.376 

0.324 

0.328 

0.329 

0.336 

0.283 

0.291 

0.294 

0.304 

0.252 

0.264 

0.269 

0.285 

0.232 

0.252 

0.261 

0.289 

0.242 

0.280 

0.299 

0.365 
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Table 31. Normalized stress intensity factors calculated at the 
midsection of an internal planar elliptic crack with 
d/h=0.25 in a plate under uniform bending M. 


V K o 


a/L Q -*-1.0 2.0 3.0 4.0 7.0 10.0 Plane 


L 0 /h 

Strain 

.05 

0.506 

0.516 

0.519 

0.520 

0.523 

0.523 

0.527 

.15 

0.535 

0. 563 

0.573 

0.578 

0.585 

0.589 

0.599 

.25 

0.5 85 

0.633 

0.652 

0.662 

0.677 

0.6 85 

0.708 

.35 

0.670 

0.751 

0.787 

0. 808 

0. 840 

0.856 

0.907 

.45 

0.895 

1.066 

1.152 

1.206 

1.298 

1.346 

1.515 


V K o 


a/L — *-1 .0 

2.0 

3.0 

V h 

V 



.05 

0.457 

0.466 

0.46 9 

.15 

0.388 

0.412 

0.421 

.25 

0.332 

0.370 

0.385 

.35 

0.285 

0.340 

0.363 

.45 

0.305 

0.332 

0.369 


4.0 7.0 10.0 Plane 

Strain 


0.471 

0.473 

0.473 

0.477 

0.425 

0.431 

0.434 

0.443 

0.393 

0.405 

0.410 

0.427 

0.377 

0.398 

0.407 

0.438 

0.395 

0.437 

0.459 

0.513 
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Table 32. Normalized stress intensity factors calculated at the 
midsection of an internal planar elliptic crack with 
d/h=0.30 in a plate under uniform bending M. 





a/L_—»-1. 0 

2.0 

3.0 

4.0 

7.0 

10.0 

Plane 

L 0 /h 








Strain 

.05 


0.600 

0.614 

0.618 

0.621 

0.624 

0.625 

0.630 

.15 


0.627 

0.666 

0.680 

0.688 

0.698 

0.703 

0.720 

.25 


0.698 

0.771 

0.801 

0.817 

0.842 

0.855 

0.895 

.35 


0. 910 

1.067 

1. 141 

1. 187 

1. 261 

1.300 

1.432 





v * 0 






a/L_— *-1. 0 

2.0 

3.0 

4.0 

7.0 

10.0 

Plane 

L 0 /h 








Strain 

.05 


0.550 

0.563 

0.567 

0.570 

0.572 

0.574 

0.579 

.15 


0.1177 

0.511 

0.524 

0.530 

0.539 

0.543 

0.557 

.25 


0.423 

0.479 

0.501 

0.514 

0.533 

0.542 

0.570 

•*35 


0.381 

0.471 

0.513 

0.538 

0.579 

0.600 

0.669 


Table 33. Normalized stress intensity factors calculated at the 
midsection of two opposite planar elliptic edge cracks 
in a plate under uniform tension N. 


a/ 

L 0 /h 

O 

t 

o 

2.0 

3.0 

4.0 

10.0 

100.0 

Plane 

Strain 

0.1 

1.060 

1.089 

1.099 

1.104 

1.113 

1.119 

1.123 

0.2 

1.009 

1.062 

1.081 

1.091 

1. 109 

1.121 

1. 125 

0.3 

0.966 

1.028 

1.065 

1.079 

1. 106 

1.124 

1.129 

0.4 

0. 929 

1.019 

1.053 

1. 072 

1. 108 

1. 131 

1. 137 

0.5 

0.902 

1.008 

1.050 

1.073 

1.118 

1.148 

1.155 

0.6 

0. 902 

1.038 

1.080 

1. 108 

1. 165 

1.203 

1.212 

0.7 

0. 929 

1.082 

1.149 

1. 186 

1.262 

1.315 

1.325 

0.8 

0.997 

1. 195 

1.284 

1.336 

1.445 

1.524 

1.539 
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Table 34. Maximum normalized stress intensity factors for 
two planar elliptic edge cracks in a plate 
under pure tension N. 


L 0 /h 

b 

a/L c 

/a — *. 0.1 

1.0 

4.0 

20.0 

Single 

Crack 

0.1 

2 

0.996 

0.986 

0.982 

0.981 

0.981 

0.1 

4 

1.072 

1.065 

1.062 

1.062 

1.062 

0.1 

10 

1.127 

1.122 

1.121 

1.120 

1.120 

0.2 

2 

0.979 

0.95 9 

0.951 

0.949 

0.949 

0.2 

4 

1 . 1 07 

1.0 90 

1.084 

1.082 

1.082 

0.2 

10 

1.221 

1.210 

1.207 

1.206 

1.206 

0.3 

2 

1.007 

0.975 

0.963 

0.961 

0.961 

0.3 

4 

1.189 

1. 160 

1.151 

1.149 

1.149 

0.3 

10 

1.382 

1.364 

1.359 

1.358 

1.358 


2 

1.051 

1.006 

0.991 

0.989 

0.989 

■ 

4 

1.2 99 

1.255 

1.243 

1.240 

1.240 


10 

1.600 

1.571 

1.563 

1.562 

1. 562 

0.5 

2 

1.106 

1.048 

1.032 

1.030 

1.030 

0.5 

4 

1.430 

1.370 

’1.354 

1.352 

1.352 

0.5 

10 

1.87 9 

1.836 

1.824 

1.822 

1.821 

0.6 

2 

1.230 

1.096 

1.0 80 

1.077 

1.077 

0.6 

4 

1.568 

1.491 

1.472 

1.469 

1.469 

0.6 

10 

2.208 

2.141 

2.124 

2.121 

2.121 

0.7 

2 

1.370 

1.120 

1.100 

1.097 

1.057 

0.7 

4 

1.694 

1.575 

1.553 

1.550 

1.550 

0>7 

10 

2.532 

2.434 

2.409 

2.405 

2.405 
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Table 35. Maximum normalized stress intensity factors for 
two planar elliptic edge cracks in a plate 
under uniform bending M. 


L 0 /h 

bi 

a/L 0 

fz —0.1 

1.0 

4.0 

20.0 

Single 

Crack 

0.1 

2 

0.874 

0.864 

0.861 

0.860 

0.860 

0.1 

4 

0.943 

0.936 

0. 934 

0.933 

0.933 

0.1 

10 

0.992 

0.988 

0. 987 

0.986 

0.986 

0.2 

2 

0.766 

0.728 

0.719 

0.718 

0.718 

0.2 

4 

0.847 

0.830 

0.825 

0.824 

0.824 

0.2 

10 

0.940 

0.930 

0.927 

0.927 

0.927 

0.3 

2 

0.751 

0.665 

0.651 

0.650 

0.650 

0.3 

4 

0. 803 

0.755 

0.745 

0.744 

0.744 

0.3 

10 

0.915 

0.900 

0.896 

0.896 

0.895 

0.4 

2 

0.792 

0.677 

0.659 

0.658 

0.656 

0.4 

4 

0.801 

0.726 

0.713 

0.711 

0.711 

0.4 

10 

0.923 

0. 902 

0.896 

0.895 

0.895 

0.5 

2 

0.826 

0.684 

0.663 

0.661 . 

0.659 

0.5 

4 

0. 834 

0.719 

0.703 

0.701 

0.700 

0.5 

10 

0.950 

0.910 

0.902 

0.901 

0.901 

0.6 

2 

0.855 

0.686 

0.662 

0.660 

0.659 

0.6 

4 

0.90 9 

0.743 

0.724 

0.722 

0.721 

0.6 

10 

0.995 

0.925 

0. 912 

0.910 

0. 910 

0.7 

2 

0.874 

0.683 

0.658 

0.655 

0.654 

0.7 

4 

0.989 

0.784 

0.761 

0.759 

0.757 

0.7 

10 

1.064 

0.956 

0.939 

0.937 

0.936 
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Table 36. Normalized stress Intensity factors calculated at the 
midsection of symmetrically located (d=0) two identical 
planar internal elliptic cracks in a plate under 
uniform tension N. 


b/a — ** 0.1 1.0 4.0 20.0 Single 


V h 

a/L 0 

Crack 

0.1 

2 

0.977 

0.976 

0.976 

0.976 

0.967 

0.1 

4 

0.987 

0.987 

0.987 

0.986 

0.987 

0.1 

10 

0.994 

0.993 

0.993 

0.993 

0.993 

0.2 

2 

0. 975 

0. 972 

0.971 

0.971 

0.971 

0.2 

4 

0.995 

0.994 

0.993 

0.993 

0.993 

0.2 

10 

1.008 

1.007 

1.007 

1.007 

1.007 

0.3 

2 

0.980 

0.980 

0.979 

0.97 9 

0.979 

0.3 

4 

1.015 

1.013 

1.012 

1.012 

1.012 

0.3 

10 

1.035 

1.034 

1.034 

1.034 

1.034 

0.4 

2 

1.007 

1.002 

1.000 

1.000 

1.000 

0.4 

4 

1.051 

1.048 

1.047 

1.047 

1.047 

0.4 

10 

1.080 

1.078 

1.078 

1.078 

1.078 

' 0.5 

2 

1.046 

1.039 

1.037 

1.036 

1.036 

0.5 

4 

1.106 

1.102 

1.101 

1.101 

1.101 

0.5 

10 

1.147 

1.145 

1.145 

1.145 

1.145 

0.6 

2 

1.108 

1.098 

1.095 

1.095 

1.095 

0.6 

4 

1.190 

1.185 

1.183 

1.183 

1.183 

0.6 

10 

1.248 

1.246 

1.245 

1.245 

1.245 

0.7 

2 

1.205 

1.192 

1.188 

1.187 

1.187 

0.7 

4 

1.321 

1.313 

1.311 

1.310 

1.310 

O'. 7 

10 

1.407 

1.404 

1.403 

1.403 

1.403 

0.8 

2 

1.367 

1.348 

1.342 

1.341 

1.341 

0.8 

4 

1.541 

1.529 

1.526 

1.525 

1.525 

0.8 

10 

1.681 

1.676 

1.674 

1.674 

1.674 

0.9 

2 

1.703 

1.672 

1.662 

1.661 

1.661 

0.9 

4 

2.007 

1.988 

1.982 

1 .981 

1.981 

0.9 

10 

2.285 

2.275 

2.272 

2.272 



72 


Table 37. Normalized stress intensity factors calculated at the 
midsection of two identical planar internal elliptic 
cracks in a plate under uniform tension N. 


V K o 

b/a — » 0. 1 1.0 4.0 20.0 


d 

L 0 /h 

.a/L 0 


0.1 

0.1 

2 

0.957 

0.821 

0.981 

0.979 

0.1 

0.1 

4 

1.002 

0.994 

0.992 

0.991 

0.1 

0.1 

10 

1.000 

0.999 

0.999 

0.999 

0.1 

0.2 

2 

1.027 

0.986 

0.977 

0.976 

0.1 

0.2 

4 

1.005 

1.002 

1.001 

1.001 

0.1 

0.2 

10 

1.018 

1.018 

1.018 

1.018 

0.1 

0.3 

2 

1.012 

0.995 

0.990 

0.989 

0.1 

0.3 

4 

1.034 

1.032 

1.031 

1.031 

0.1 

0.3 

10 





0.1 

0.35 

2 

1.021 

1.010 

1.005 

1.004 

0.1 

0.35 

4 

1.058 

1.056 

1.055 

1.055 

0.1 

0.35 

10 





0.2 

0.1 

2 

0.947 

0.783 

m 


0.2 

0.1 

4 

1.006 

0.995 



0.2 

0.1 

10 

1.003 

1.002 

mm 


0.2 

0.2 

2 

1.052 

0.993 

0.981 

0.979 

0.2 

0.2 

4 

1.020 

1.016 

1.014 

1.014 

0.2 

0.2 

10 

1.039 

1.038 

1.038 

1.038 

- 0.2 

0.25 

2 

1.056 

1.003 

0.994 

0.992 

0.2 

0.25 

4 

1.044 

1.040 

1.040 

1.039 

0.2 

0.25 

10 

1.073 

1.073 

1.073 

1.073 

0.3 

0.1 

2 

0.933 

0.760 

0.980 

0.975 

0.3 

0.1 

4 

1.01 8 

1.003 

0.999 

0.998 

0.3 

0.1 

10 

1.014 

1.013 

1.013 

1.013 

0.3 

KEI 

2 

0.814 

0.985 

0.988 

0.985 

0.3 

maS 

4 

1.033 

1.024 

1.022 

1.021 

0.3 

KSS 

10 

1.047 

1.046 

1.046 

1.046 
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Table 37 - Cont 


V K o 


b/a _*» 0. 1 1.0 4.0 20 .0 


d 

L c /h 

a / L 0 


0.1 

0.1 

2 

0.95 9 

0.830 

0.982 

0.979 

0.1 

0.1 

4 

1.002 

0.993 

0.991 

0.991 

0.1 

0.1 

10 

0.999 

0.999 

0.999 

0.999 

0.1 

0.2 

2 

1.020 

0.984 

0.978 

0.976 

0.1 

0.2 

4 

1.002 

1.000 

0.999 

0.999 

0.1 

0.2 

10 

1.014 

1.014 

1.014 

1.014 

0.1 

0.3 

2 

1.004 

0.989 

0.985 

0.984 

0.1 

0.3 

4 

1.021 

1.020 

1.019 

1.019 

0.1 

0.3 

10 





0.1 

0.35 

2 

1.006 

0.997 

0.993 

0.993 

0.1 

0.35 

4 

1.036 

1.035 

1.034 

1.034 

0.1 

0.35 

10 





0.2 

0.1 

2 

0.950 

0.799 

0.979 

0.976 


0.1 

4 

1.004 

0.994 

0.992 

0.991 

0.2 

0. 1 

10 

1.001 

1.001 

1.001 

1.001 

mm 

0.2 

2 

1.038 

0.987 

0.977 

0.976 

0.2 

0.2 

4 

1.010 

1.007 

1.006 

1.005 


0.2 

10 

1.026 

1.026 

1.026 

1.026 


0.25 

2 

1.034 

0.990 

0.983 

0.982 

SB 

0.25 

4 

1.023 

1.021 

1.020 

1.019 

0.2 

0.25 

10 

1.046 

1.046 

1.046 

1.0 46 

0.3 

0.1 

2 

0.935 

0.778 

0.978 

0.974 

0.3 

0.1 

4 

1.013 

0.999 

0.996 

0.995 

0.3 

0.1 

10 

1.009 

1.009 

1.008 

1.008 

0.3 

0.15 

2 

0 . 829 

0.977 

0. 579 

0.977 

0.3 

0.15 

4 

1.01 9 

1.011 

1.008 

1.008 

0.3 

0.15 

10 

1.030 

1.030 

1.030 

1.030 
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Table 38. Normalized stress intensity factors calculated at 
the midsection of three identical planar internal 
elliptic cracks in a plate under uniform tension N, 


K A /Ko for the middle crack 


d 

L 0 /h 

a/L 0 

• 

o 

1 

nJ 

\ 

.0 

1.0 

4.0 

20.0 

0.1 

0.1 

2 

0.983 

0.981 

0.979 

0.979 

0.1 

0.1 

4 

0.994 

0.992 

0.992 

0.991 

0.1 

0.2 

2 

0.985 

0.980 

0.977 

0.976 

0.1 

0.2 

4 

1.007 

1.004 

1.002 

1.001 

0.1 

0.3 

2 

1.006 

0.997 

0.991 

0.989 

0.1 

0.3 

4 

1.040 

1.035 

1.032 

1.031 

0.2 

0.1 

2 

0.981 

0.978 

0.976 

0.975 

0.2 

0.1 

4 

0.995 

0.994 

0.992 

0.992 

0.2 

0.2 

2 

.0. 993 

0.986 

0.981 

0. 979 

0.2 

0.2 

4 

1.023 

1.018 

1.015 

1.014 

0.3 

0.1 

2 

0.983 

0.979 

0.976 

0.975 

0.3 

0.1 

4 

1.003 

1.001 

0.999 

0.998 


K b /K q for 

the middle crack 



0.1 

0.1 

2 1 

0.984 

0.981 

0.980 

0.979 

0.1 

0.1 

4 

0.993 

0.992 

0.991 

0.991 

0.1 

0.2 

2 

0.984 

0.980 

0.977 

0.576 

0.1 

0.2 

4 

1.003 

1.001 

0.999 

0.999 

0.1 

0.3 

2 

0.997 

0.990 

0.986 

0.984 

0.1 

0.3 

4 

1.026 

1.022 

1.020 

1.019 

0.2 

0.1 

2 

0.981 

0.578 

0.977 

0.576 

0.2 

0.1 

4 

0.954 

0.993 

0.992 

0.991 

0.2 

0.2 

2 

0.987 

0.981 

0.977 

0.976 

0.2 

0.2 

4 

1.012 

1 .009 

1.006 

1.005 

0.3 

0.1 

2 

0.981 

0.577 

0.975 

0.974 

0.3 

0.1 

4 

0.999 

0.997 

0.996 

0.995 
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Table 38 -Cont 


K a / K q for outer crack 


d 

L 0 /h 

a / L 0 

b/a — •> 0.1 

1.0 

4.0 

20.0 

0.1 

0.1 

2 

0.982 

0.980 

0.979 

0.979 

0.1 

0.1 

4 

0.993 

0.992 

0.991 

0.991 

0.1 

0.2 

2 

0.982 

0.978 

0.976 

0.976 

0.1 

0.2 

4 

1.005 

1.003 

1.002 

1.001 

0.1 

0.3 

2 

0.999 

0.994 

0.991 

0.989 

0.1 

0.3 

4 

1.037 

1.034 

1.032 

1.031 

0.2 

0.1 

2 

0.979 

0 . 977 . 

0.976 

0.975 

0.2 

0.1 

4 

0.994 

0.993 

0.992 

0.992 

0.2 

0.2 

2 

0.988 

0.983 

0.980 

0.979 

0.2 

0.2 

4 

1.019 

1.016 

1.014 

1.014 

0.3 

0.1 

2 

0.980 

0.978 

0.976 

0.975 

0.3 

0.1 

4 

1.001 

1.000 

0.999 

0.998 


K s / K 0 for outer crack 




0.1 

0.1 

2 

0.982 

0.980 

0.980 

0.979 

0.1 

0.1 

4 

0.993 

0.992 

0.991 

0.991 

0.1 

0.2 

2 

0.981 

0.978 

0.976 

0.976 

0.1 

0.2 

4 

1.002 

1.000 

0.999 

0.999 

0.1 

0.3 

2 

0.992 

0.988 

0.985 

0.984 

0.1 

0.3 

4 

1.023 

1.021 

1.020 

1.019 

0.2 

0.1 

2 

0.979 

0.978 

0.976 

0.976 

0.2 

0.1 

4 

0.993 

0.992 

0.991 

0.991 

0.2 

0.2 

2 

0.983 

0.979 

0.977 

0.976 

0.2 

0.2 

4 

1.010 

1.007 

1.006 

1.005 

0.3 

0.1 

2 

0.97 8 

0.976 

0.974 

0.974 

0.3 

0.1 

4 

0.998 

0.996 

0.995 

0.995 
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Table 39. Normalized stress intensity factors on the crack front 
for an internal elliptic crack with d/hrO.20, L /h=0.45, 
a/L 0 =4. 


X 

<Wh 

(k b /k 0 

0.99 

0.426 

0.411 

0.90 

0.723 

0.707 

0.80 

0.865 

0.831 

0.70 

0. 974 

0. 916 

0.60 

1.065 

0.980 

0.50 

1.143 

1.028 

0.40 

1.211 

1.06 8 

0.30 

1.266 

1.096 

0.20 

1.307 

. 1.116 

0.10 

1.333 

1.128 

0.00 

1.342 

1.132 


N (K B /K o^M (K B /K o^M 


0.186 

0.151 

0.357 

0.218 

0.454 

0.229 

0.531 

0.233 

0.595 

0.233 

0.64 9 

0.233 

0.697 

0.233 

0.735 

0.233 

0.763 

- 0.233 

0.781 

0.233 

0.787 

0.232 


77 



Table 40. Normalized stress intensity factors on the crack front 
for an internal elliptic crack with d/h=0.10, L o /h=0.45, 
a/L 0 =4. 


X 

<Wn 

<Wn 

<Wm 

< K B /K o>M 

0.99 

0.411 

0.393 

0.097 

0.065 

0.90 

0.692 

0.686 

0.205 

0. 072 

0.80 

0.816 

0.80 4 

0.269 

0.057 

0.70 

0. 902 

0.883 

0.318 

0.041 

0.60 

0. 967 

0.941 

0.357 

0.028 

0.50 

1.018 

0.984 

0.388 

0.016 

0.40 

1.061 

1.019 

0.414 

0.007 

0.30 

1.092 

1.044 

0.433 

0.000 

0.20 

1.114 

1.061 

0.446 

-0.005 

0.10 

1.128 

1.072 

0.455 

-0.008 

0.00 

1.132 

1.075 

0.457 

-0.009 
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Table 41. Normalized stress intensity factors on the crack front 
for an elliptic edge crack with Lg/hsO.45, a/L 0 =4. 


X 

<w 

<w 

0.99 

0.762 

0.646 

0.90 

0.920 

0.696 

0. 80 

1.027 

0.711 

0.70 

1.084 

0.701 

0.60 

1.130 

0.695 

0.50 

1.166 

0.689 

0.40 

1.197 

0.686 

0.30 

1.214 

0.680 

0.20 

1.227 

0.676 

0.10 

1.237 

0.675 

0.00 

1.240 

0.675 
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FIGURES 
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Figure 3. 


(b) 


Geometry and notation for an internal crack. 
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I 


Figure 5. Geometry and notation for two identical cracks. 



Figure 6. Geometry and notation for three identical cracks. 
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1: 

L /hsO.65, 

a/L 0 =7 

2: 

L /h=0. 65 ) 

a/L o=3 

3: 

L /hs0.55, 

a/L =7 


L°/h=0. 55 , 

a/L =3 

5 : 

L°/h=0. 45 , 

a/L o=7 

6 : 

L /h=0.J»5, 

a/L°=3 

7: 

L°/h=0.35, 

a/L =7 

8: 

L /h=0. 35, 

a/L =3 

9: 

L°/h=0.25, 

a/L o=3 


Figure 7. Effect of d/h ratio on the normalized stress intensity 
factors calculated at the midsection of various 
internal elliptic cracks under pure tension. 



0.05 0.1 


015 


0-2 


0.25 03 d/h 


1: 

L /h=0. 35, 

a/L =3, 

2: 

L o /h=0.35, 

a/L 0 =7, 

3: 

L /h=0.65, 

a/L =3, 

4: 

L /h=0.55, 

a/L =3, 

5: 

L /h=0.45, 

a/L =3, 

6: 

L /h=0.45, 

a/L =7, 

7: 

L /h=0.55, 

a/L°=7 , 

8: 

L /h=0.65, 

a/L =7, 

9: 

L o /h=0.25, 

a/L = 3 . 


Figure 8. Effect of d/h ratio on the normalized stress intensity 
factors calculated at the midsection of various 
internal elliptic cracks under pure tension. 



0-05 0.1 0.15 0 2 0 25 


1 

L /h=0.25, 

a/L 0 =3 

2 

L /h=0.25, 

a/L Q =7 

3 

L /h=0.35, 

a/L 0 =3 

4 

L /h=0. 35, 

a/L 0 =7 

5 

L /h=0.45, 

a/L =3 

6 

L /h=0.45, 

a/L 0 =7 

7 

L /h=0.55, 

a/L 0 =3 

8 

L /h=0.55, 

a/L 0 =7 

9 

L /h=0.65, 

a/L 0 =3 

10 

L_/h=0.65, 

a/L 0 =7 

11 

L_/h=0. 15, 

a/L =3 

12 

L /h=0.15, 

a/L 0 =7 


0.3 cl/ h 


Figure 9 


Effect of d/h ratio on the normalized stress intensity 
factors calculated at the midsection of various 
internal elliptic cracks under pure bending. 




0 05 o.l 


0.15 


0.2 0.25 0-3 d/h 


Is L o /h=0.15, a/L =3, 
2: L°/h=0.25, a/IT=3, 

3: L /h=0.35, a/L =3, 

4: LVh=0.45, a/L =3, 

5: L /h=0. 55 , a/L =3, 

6: L /h=0.65, a/L =3, 


Figure 10. Effect of d/h ratio on the normalized stress intensity 
factors calculated at the midsection of various 
internal elliptic cracks under pure bending. 





0.05 



Figure 11. 


Effect of L Q /h ratio on the normalized stress 
intensity factors calculated at the midsection of 
internal elliptic cracks with a/L 0 =4 and 
different d/h ratios under pure tension. 
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Figure 16. Effect of a/L 0 ratio on the normalized stress 

intensity factors calculated at the midsection of 
internal elliptic cracks with d/h=0.15 and 
different L_/h ratios under pure bending. 


















K/Ko 


1 8 - 


1.6 - 


1 4 ■ 

j 

i 

1.2 - 
1.0 - 


a/L o =10 



01 0.3 0-5 0.7 1 u/h 


— Line spring 

Ref. [19] 

Figure 20. Comparison of the normalized stress intensity 
factors calculated in this study at x=0 for an 
internal planar elliptic crack m a plate under 
uniform tension and with the corresponding 
values given in Ref. [19]. 


99 




S3 




— K/hL 



Figure 22. Normalized stress intensity factors calculated at 
the midsection of two opposite planar elliptic 
edge cracks ih a plate under uniform tension. 



Figure 23. Comparison of the stress intensity factors calculated 
in this study at the midsection of a symmetrically 
located internal(d/h=0) planar elliptic crack in a 
plate under uniform tension with the corresponding 
results given in [20]. L Q /h=. 57, a/L Q =1.25, x=cos $. 
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APPENDIX I 


The dimensionless quantities used in the formulation 
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APPENDIX II 


Procedure for finding coefficients of the shape functions (g t 
and g^) are as follows. Loading should be seperated as pure tension 
and pure bending. For pure tension case 


m = 0 . 


By using (8), (9) and (10a) 


K(s) = \fh *g t (s) , 

K(s) = \fh9 'J*3 Jb 1 (s) 2(i " 1) , 

iml 

K(s) = flr V^LCxTT Z b i (s) 2 ^ 1 “ 1 ^ , 

Csl 

=Zb i( s) 2 ^-i) 

C y 1TL(x^ ) Ui 


Values of the left handside can easily be found for a wide range of 
s values in related literature, which enables us to create n 
equations with b^s as the n unknowns. Solution of this system is 
the coefficient of corresponding shape function. 


Same procedure is followed under pure bending to find c^'s 
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APPENDIX III 


Solution of the Integral Equations 

i- Gauss-Chebyshev Closed Type Integration Formula 
The solution of the integral equation which only has the Cauchy 
type singularity 
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-1 < x < 1 
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and the additional condition 


J $(t) dt = A , 

-I 


becomes 
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One may find the detailed derivation of the formula in [4]. 


When there is a system of integral equations with more than one 
unknown functions to be solved, the above procedure can be used by 
applying the formula to each integration seperately. For example 
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-1 < X < 1 
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J ♦^t) dt = A 
-I 
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♦ 2 (t) ^ ~ ^2 


> 


can be solved as 


♦ t (t) = (1 - t 1 2 )“ 1/2 gl (t) , 

♦ 2 (s) = (1 - s 2 r 1/2 g 2 (s) , 


s has been Introduced to aefine different collocation points for 
each functions. 
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ri-i *| 
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Above solution can also be used for any number of unknown functions 
in the same manner. 


ii- Order of Singularity in the Integral Equation 

At. first glance, it seems that (31) has singularities other 
than Cauchy type. But if we take a closer look to the equation, we 
will see that modified Bessel function of the second kind has the 
following property for relatively small values of the argument. 


K r (x)~ a^Cr-Oix" 1, , 

K 2 (x) 2 x - ^ 

By using above property, we can say that 
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which gives 


11m 

(t-x)— 0 
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This proves that (31) has only Cauchy type singularity and can be 
solved by Gauss-Chebyshev closed type integration formula. 


Integrals on the interval (-1 , x) can be written on (-1 , +1) 
by redifining the kernels by the use of a step function H(y) as 
follows: 


x 



H(x-t) f(t) dt 


where, 


H(y) 


0 if y < 0 

1 if y >, 0 . 
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APPENDIX IV 


Numerical Integration 

Gauss-Chebyshev closed type integration formula is used for 
numerical integration purposes. Any integral in the interval 
(-1 , +1) can be expressed as, 

» r 1 V 1 1 

— f(x«) + £ fCx,) +— f(x n ) 

n-1 L 2 Ua. 2 J 


I 

f(x) dx 
'/l -x 2 


If an integration in the interval (-1 , x) is desired, above 
formula can again be used by multiplying the integrand with the H 
function defined in Appendix III. 

Number of points (n) in the integration should be much higher 
than the number of points used in the solution of integral equations 
if one desires to evaluate the stress intensity factor or the crack 
front at any point accurately. Values of the functions to be 
integrated at any point can be found by interpolating the function 
between the two surrounding known values found from the solution of 
integral equations. 
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